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because this oxidizer's stored aboard 


Nitrogen Tetroxide provides immediate, hypergolic igni- 
tion with amine fuels. Even after 9 years of storage, 
nitrogen tetroxide is ready for instant use. Combustion 
efficiencies closely approach the theoretical 99%. 
N2O,—Requires no refrigeration—can be stored indefi- 
nitely in missiles at launching sites—Can be used with 
most fuels—including those containing carbon.—Elimi- 
nates rough starts—fast reaction prevents accumulation 


llied 
hemical 


of propellants in thrust chambers.—Allows throttleable 
control of motors. 

We'll gladly supply technical literature, including: a 
59-page Product Bulletin, and a brochure entitled “Large 
Scale Handling of Nitrogen Tetroxide.” 

For specifications and local offices, see our insert in 
Chemical Materials Catalog, pages 475-482 and in Chemi- 
cal Week Buyers Guide, pages 37-44. 
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will The Mission 

Succeed? 


(Dynamic Accuracy Test 


racy Test System) tells the interceptor com- erly operating weapon control systems are 


Flight-line checkout by DATS (Dynamic Accu- It makes certain that only aircraft with prop- 7 
mander whether his aircraft and weapon control sent on missions. Based on a building-block __ . ; 


systems are completely ready for a successful design employing the highest reliability fac- 

mission. As a result of field evaluation tests, tors, a mechanical programming device and 

showing the effectiveness of DATS in improv- self-test capability, DATS utilizes a series of 
ing weapon control performance, RCA has synthesized attack runs typical of mission ie 
been awarded an Air Force production con- conditions. DATS could be made applicable 


tract. Developed by RCA’s Airborne Systems A. 
Division, Defense Electronic Products, : 
Camden, New Jersey, DATS is a new ap- © 
proach to the evaluation of — readiness. 


to many interceptor types of aircraft. 
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mechanics, high speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
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to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
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The national security need for Space Technology Leadership 
Space is a medium in which many military missions can be most effectively performed @ The U. S. arsenal of ballistic missiles—in 


being and forthcoming — the Air Force Thor, Atlas, Titan, and Minuteman, the Navy Polaris, and the Army Jupiter, are all designed 
to deter the outbreak of a nuclear World War II! or to retaliate overwhelmingly if it should occur @ If our ballistic missiles are to 
realize their greatest potential in carrying out their dual task, they must be supported by a number of companion space systems 
for such missions as early warning, reconnaissance, communications, navigation, weather forecasting. Space Technology Labora- 
tories is proud of its contributions to the national space effort as a principal contractor in carrying out major programs for the Air 
Force Ballistic Missile Division, National Aeronautics and Space Administration, and Advanced Research Projects Agency @ The 
increasing scope of STL’s activity is opening up exceptional opportunities for the exceptional scientist and engineer, who will find 
creativity given encouragement and recognition in an organization synonymous with Space Technology Leadership. Resumes and 


inquiries will receive meticulous attention. 


SPACE TECHNOLOGY LABORATORIES, INC. p.o. 80x 95005c, Los ANGELES 45, CALIFORNIA 


a subsidiary of Thompson Ramo Wooldridge Inc. 


Cape Canaveral * Manchester, England * Singapore * Hawaii 
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Is the “S” curve our 
best missile defense? 


There have been many theories advanced on techniques 
of ICBM intercept, including those that postulate the 
use of an anti-missile missile. Some studies favor a head- 
on course for the weapon during terminal phase because 
any angle other than 180° from warhead flight allows too 
much error-margin, and 360° allows too much time lag 
for the “chase.” 


@ To give the anti-missile as much time as possible to seek the warhead, 
and lower the maneuverability requirements, it should be in the air ready 
to intercept while computers translate the acquisition signals into a plotted 
line of approach. Thus the anti-weapon can be closing the altitude gap 
before it “knows” the exact approach vector. The path described by this 


anti-missile would take the form of an “S” curve. os - 


@ Yet, there may be other, more effective means of intercept. At Loral, 
where fresh attacks on formalized thinking result in new concepts for 
defense, the broad approach to problems such as this is backed in depth 
by advanced research in such fields as electronic scanning and electro- 
optics. And from these probings come designs that our production teams 
turn into operational defense hardware. 


@ If you are a senior engineer or scientist interested in being part of this — 
rapidly-expanding thinking group, write LORAL Electronics Corporation, — 
New York 72, New York. 
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publish selected papers offered at the specialized 
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the same subject. Beginning with Volume 2, each 

~ volume is expected to be published within four months 
poe _ after the papers are received. 


px: _ At this time, when the quantity of scientific papers is 
ca rapidly exceeding the capacities of the customary media 
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Structures Considerations in Design 
for Space Boosters 


PAUL E. SANDORFF'! 


Massachusetts Institute of Technology 
| Cambridge, Mass. 


Structures aspects of the design of space vehicles are studied from the viewpoint of how they in- 


oe ~~ - fluence and are influenced by other technical problems of the space mission. Booster size and con- 


- figuration are given special attention. 
larger vehicles than those currently contemplated, and that a systems approach to the space effort 
“es re, calls for strong emphasis on the structures contribution. Sie a 


a, 


HE OBJECT of this paper is to describe the ways in 
which structures problems influence and are influenced 
by the other aspects of early space vehicle design. 

To achieve this aim, the nature of the major loading con- 
ditions and how they may vary with vehicle size, configuration 
and design parameters is first discussed. Structural design 
methods which are available to support the major vehicle 
loads are next compared, using structural indexes to assess 
future design problems. Functional components of the 
structure are analyzed to determine how weight may vary 
with vehicle design parameters. Dead weight data available 
on existent designs are summarized, and the trends noted. 
Finally, in support of the predictions derived analytically, 
a family of large booster vehicle designs is presented. 

Although at this time all space vehicles are chemical 
rockets which take off from the surface of Earth, the modular 
concept of orbital assembly, pilotless rendezvous techniques 
and orbital refueling may soon provide preliminary steps 
toward orbital erection and satellite-based, high performance 
interplanetary vehicles. Such measures are only another 
form of staging, where the first stage is comprised of multiple 
ascents of a device of limited size instead of a single large 
booster, and as such they deserve re-evaluation. A particular 
objective of this paper therefore is to predict the effect of size 
increase on the structures aspects of large liquid propellant 
booster rockets. 


System Design and the Structural Requirements 


The interrelationship between vehicle design considerations 
and major structural loading conditions is considered here 
from the point of view - minimising structural problems 


through design, 


Handling requirements may be quite significant to the 
structural design, and problems of transportation from fac- 
tory to launch pad have strongly influenced contemporary 
booster dimensions. A serious venture into space demands 
re-evaluation of our Earth-side methods and facilities from a 
systems point of view. 

Loads which occur on the launching pad include a ground 


Presented at the ARS 13th Annual Meeting, New York, N. Y., 
Noy. 17-21, 1958. Revised, August 1960. 

| Associate Professor, Department of Aeronautics and Astro- 
nautics. Member ARS. 
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It is concluded that structures considerations favor much 


wind condition, important when the vehicle is supported only 
at the base. In addition to the static loading, von Karman 
vortex excitation can induce the low frequency cantilever 
bending modes, possibly leading to a fatigue condition (1).? 
Consideration should be given to structural design with ex- 
ternal stiffeners which serve as spoilers. The loading con- 
dition may not be a source of great concern for vehicles 
of smaller fineness ratio than contemporary designs. 


Thrust Loading 


A major portion of the structural weight results from the 
design condition owing to thrust, which is reacted by inertial 
loading throughout the vehicle. Maximum acceleration for 
any stage occurs just prior to burnout and establishes design 
values for thrust transmitted to forward structure and pay- 
load. 

In propellant tanks, hydrostatic pressure is the mechanism 
whereby the fluid is accelerated, and the thrust condition 
therefore establishes tankage design pressures. In the first 
stage this pressure component is maximum at the instant of 
takeoff; for upper stages the hydrostatic loads, like the 
other thrust loads, will usually be established at booster burn- 
out. Hydrostatic pressure requirements may be raised 
significantly by dynamic effects arising from longitudinal 
thrust variations and the symmetrical surging associated with 
fluid compressibility and tank wall elasticity. An overload 
factor approaching 2 is possible, especially if the time for 
acceleration to build up is less than the period of the surging 
mode. The undamped natural period for this dynamic sys- 
tem in a simple cylindrical tank is approximately (2) 


For contemporary applications, this period is already appreci- 
able (0.2 sec for the Saturn booster, for example) compared 
to the time for thrust buildup at takeoff. In addition to the 
simple surging overload is the possibility of self-supported 
dynamic response owing to coupling of pump inlet pressure 
with engine thrust. In liquid hydrogen tankage, the high 
fluid compressibility will make surging effects slower and 
therefore more difficult to control. 

Little can be done to reduce the thrust requirements for 
the booster stage of Earth takeoff vehicles. In fact, slightly 


2 Numbers in parentheses indicate References at end of paper. 
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higher thrust factors might be anticipated as booster size is 
increased, because larger vehicles generally have larger 
ballistic coefficients, and drag becomes less important in 
determining the overall performance optimum. Upper stage 
thrust loading, however, may be strongly influenced by over- 
all design considerations: 

1 Lower burnout acceleration for the previous stage 
resulting from improvement in specific impulse or the use of 
multiple staging. 

2 The use of low thrust, long duration rocket engines as a 

better overall compromise with the trajectory penalties due 
to gravitational losses. 
Losses associated with the optimum trajectory diminish in 
magnitude as the vehicle speed approaches and exceeds 
orbital velocity. For a liquid propellant rocket operating 
from orbit, for example, engine mass considerations alone 
result in a thrust level of perhaps 4 G as the best compromise 
between the mass ratio penalty of a larger heavier engine and 
the gravitational loss of a lower thrust, longer duration 
engine. Vehicles which are to be erected or filled in orbit will 
strike an even lower compromise thrust level, because this 
thrust condition will also establish structural mass of the 
vehicle. 


Pressurization 


In addition to hydrostatic pressures already mentioned, 
internal pressure loads may be expected for the following 
reasons: 

1 In tankage structure, for propellant handling purposes. 
In membrane structure, for stabilization purposes. 
3 In eabins of manned vehicles, to provide artificial en- 
vironment. 

4 In cargo and instrument compartments, for functional 
purposes, such as temperature control. 

5 In propulsion system components: Combustion cham- 
bers, heat exchangers, piping, etc. 

6 In compressed gas storage vessels. 

Clearly, space vehicle structure will be largely pressure 
vessel design. 


Maneuvering and control conditions provide the only ex- 
ternal transverse loads of any magnitude. For space opera- 
tions the time scale can be of such an extent that maneuvering 
requirements per se will establish only minor loads. Control 
loads may bear further examination, depending on the func- 
tions of the control system, detail problems of control which 
are not yet clearly defined and the way in which the science 
may develop to meet them. 

One such problem concerns the interaction of propellant 
sloshing with engine thrust vectoring, vehicle bending and 
control dynamics (3,4). The natural frequency of the side- 
to-side sloshing of fluid in a circular cylindrical tank with the 
free surface more than one diameter above the bottom may 
be roughly predicted by the formula : 


Sa viosh = 1.7V/n/D 


where 


n = acceleration factor, Earth G 
D = tank diameter, ft 


This is an empirical relation based on small scale tests. It 
indicates extremely low frequencies for the large low accelera- 
tion vehicles which might be visualized for interplanetary 
journeys. The forces involved depend on the mass of the 
sloshing fluid, and for geometrically similar amplitudes this 
varies as the cube of the tank diameter (3). Consequently, 
large size tankage will be possible only if sloshing can be con- 
trolled. 

This problem may be attacked in several ways; in any 


case, consideration of all dynamical effects simultaneously, 
analytically or experimentally, is long overdue (5,6). Nelson 
(5) has recently demonstrated that the simple pendulum 
analogy generally used to represent propellant sloshing may 
be inadequate; sloshing raises the natural frequency and 
effective damping of structural modes, and drainage of fluid 
from the tank produces an important negative damping effcct 


on the sloshing phenomenon itself. His work was also 
directed to provide generalized relations for pressure <is- 
tribution and system response, on the basis of which presstire 
pickups and a servo system might be designed to introduce 


artificial damping of sloshing by out-of-phase thrust vecior 
control. This type of artificial damping is commonly used 


to reduce overall vehicle bending disturbances, and was 


first applied in the Vanguard (7). 
Control of sloshing would be much easier if the freque: «y 


of the vehicle autopilot system could be made low compare to 
the frequency of the lowest sloshing mode. This can cr- 
tainly be done while still meeting the trajectory limitations 
and maneuvering requirements of spaceflight in designs fo :c- 
seen for chemical rockets. For exotic propulsion schetiics 


which result in very low vehicle acceleration, the interaction 
of fluid stores with structure and vehicle response may ag:in 
exert important and as yet unexplored limitations on systim 


Booster Flight Loads 


An important design load problem for Earth takeoff 
rockets is the complex one involving aerodynamic disturbaiice 
owing to lateral winds acting on the ascending vehicle and the 


dynamic response of the vehicle and its control system. A 
typical design condition stipulates steadily increasing lateral 


wind (wind shear condition) tending to 150 fps at 40,000-ft 


altitude, then diminishing; superimposed on this profile is a 


triangularly shaped spike (jetstream) of an additional 150-{ps 


strength, occupying perhaps 5000 ft of vertical distance in 


the 30,000- to 40,000-ft altitude zone. It is in this altitude 
range that maximum dynamic pressure and maximum 4q/? 


occur for large liquid propellant rocket boosters; a typical 
mission for such devices would result in a velocity of perhaps 
_ 2000 fps at an inclination 30 to 40 deg above the horizontal at 
40,000-ft altitude. 
_ Of primary significance in determining the resultant struc- 
tural loads is the guidance philosophy. Control of the 
vehicle flight path may be based on the attempt to hold a 
constant or programmed attitude with respect to inertial 
space; even with perfect control, this results in substantial 
aerodynamic forces. Lower loads and a better trajectory 
would be obtained with a system which attempts to hold, in- 
stead of a specified attitude, a specified flight path. The 
minimum aerodynamic load is obtained by exercising control 
with respect to angle of attack (weathercock stability), with, 
however, considerable trajectory displacement. Geissler (4), 
in an unusually thorough treatment of this problem, has 
shown how by obtaining angle of attack information in addi- 
tion to inertial attitude and by using programmed gain fac- 
tors, a swivel engine system which is maneuverable with 
respect to the disturbance may be designed to achieve charac- 
teristics varying from “minimum drift control” to “minimum 
load control,’ the latter affording significant reduction in 
structural and control requirements. 

Equal in importance in determining structural loads is the 
capability of the control system, for if it is marginal in capac- 
ity or in speed of response, high structural loads and even di- 
vergence may result. But practical problems limit the control 
system design: Favorable aerodynamic characteristics usu- 
ally entail compromises and weight penalties, as does actuator 
capability. Control system dynamical response is limited by 
structural vibratory modes, which the gyro and accelerometer 
sensors confuse with rigid body motion. To avoid coupling 
effects which may result in saturation or instability, the lower 
modes may be damped by phase shifting of thrust response 
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(7), high frequency signals blocked with filtering networks, 
and the vehicle autopilot response consequently held to 
perhaps a fifth of the fundamental body bending frequency. 
There is little doubt, however, that this ratio can be improved, 
with greater sophistication in control system design, for 
example, by employing multiple sensors or strain gage net- 
works which permit identification and rejection of the effects 
arising from structural deformation. 

How are the control system requirements affected by size 
and shape of the booster? For large boosters these require- 
ments are established chiefly by the jetstream, which the 
vehicle following a reasonable trajectory passes through in 
perliaps 5 see regardless of size. Assuming vehicle pitch 
ang!e to be kept small compared to the angle of attack, the 
moment impulse developed by the aerodynamic dis- 
tur)ance is 


Jf = f Cng (Fp: (kL)dr 


Fo: cone bodies the coefficient Cm, Varies slightly with L/D, 
giving some incentive for lower fineness ratio, that will be 
disregarded for the time being. Adopting Cm, as invariant 
on the basis indicated and using k to indicate proportionality, 
th aerodynamic input over any selected time increment 
varies among similar vehicles as - 


f = f vdr 
v 


The rate of response required of similarly effective control 
systems will be determined by the rate of pitch which the 
vehicle develops, for if it pitches at a slow rate, the corrective 
maneuver may also be accomplished leisurely—assuming 
always that the total excursion will not be permitted to be 
large. Making an approximation for the instantaneous 
vehicle moment of inertia 


Iver = kW pl? = kppD*L 


the minimum required natural frequency of the vehicle auto- 
pilot system is seen to vary as 


ontro min ppl 
oh 


where the quantity JS (u/v)v%dr has been assumed invariant 
for vehicles following similar trajectories and exposed to 
similar aerodynamic disturbances. In large vehicles, there- 
fore, the control system may operate much more slowly and 
still maintain control. 

The problem of interaction of the control system with 
structural modes depends on whether vehicle bending fre- 
quencies are similarly reduced. For clustered tank designs, 
the bending rigidity is unusually low, and a worsening of the 
situation can be avoided only by accepting sizable weight 
penalties. Assuming that single cylindrical shell tankage is 
to be used, the overall structural bending frequency may be 
assumed to vary roughly as 


bend. > kV El 
If the structural weight fraction for the shell is to remain the 
sume among the vehicles considered, that is 
W. _ kp.DLt 
the shell moment of inertia may be estimated as 
Wp D? 


= constant 


Substituting back in the frequency relation with the material 
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properties considered constant 

fren bena. = kD/L? 
In comparison with the control frequency relation, it must be 
expected that increase in vehicle diameter will relieve the 
problem of interaction betweer structural modes and auto- 
pilot. Such an effect appears inescapable as the gross size of 
vehicles is increased. Reduction in propellant density, how- 
ever, shifts the required control frequency up into the range 
of frequencies of structural modes. 

The magnitude of the control deflection necessary to avoid 
exceeding any specified pitching response is established by the 
requirement that the system provide an impulse sufficient to 
overcome the aerodynamic input. The available control 
impulse will vary approximately as 


W 
f = Jf sar = kn; = prD*L Bdr 
Wp 


Kquating to Mecrolt and solving 
1 Wpl = 


ni Wi pp 


In the wind-shear type of disturbance, the control response 
will keep up with the aerodynamic input, showing a lag which 
could be used as a measure of effectiveness. Among equally 
effective systems the required control deflection is therefore 
determined by the instantaneous aerodynamic angle of 
attack w/v, and since this will vary in time about the same for 
all vehicles, timewise integration of the preceding relation 
leads to time factors on both sides which are proportionate to 
each other. Then 


For indicial disturbances, however, the aerodynamic im- 
pulse integrates over a fixed time scale, whereas the control 
impulse integrates over a period which was found to vary as 
ppL*. In this case the maximum required control deflection 
will vary for different vehicle sizes and shapes approximately 


as 


1 Wel 
k ne Wi pp _ Wr 
max req. ppL? nN; Wi pp*L? 


Since the control system period increases with L?, it is con- 
ceivable that for very large vehicles even the jetstream dis- 
turbance may be considered an indicial input. 

In summary, large vehicles will be sluggish in motion, but 
it appears that the control system problem will actually be 
alleviated by increase in vehicle size. This result is in 
accord with an earlier demonstration by Hebeler (8) based on 
a study of related designs of size ranging from 105 to 10° lbm 
gross, which were “flown” through a synthetic design wind 
profile on a computer to obtain the maximum pitch angle vs. 
gross weight; the conclusion reached was that very large 
booster vehicles will go through the most severe aerodynamic 
disturbance “like the Queen Mary.” 

Propellant sloshing, which has been discussed in the fore- 
going, is of considerable concern in contemporary boosters, 
because sloshing frequencies vary over a range which ap- 
proaches and may bracket the value that other considerations 
dictate for the vehicle autopilot system (3,4). However, 
sloshing frequencies vary approximately as 1/+/D, whereas, 
as has been seen, the required vehicle autopilot natural fre- 
quency varies approximately as 1/L?. Consequently the 
problem of dynamic interaction is expected to diminish as 
the two frequencies diverge in very large boosters. Active 
damping measures may then serve to control incidental 
sloshing, such as properly phased thrust vector control as 
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previously suggested, or some other wave or circulation- 
producing system. 

This reasoning is based on a study of propellants of density 
similar to those in use today. The effect of propellant density 
p, enters in all the relations derived here in such a way as to 
make the problems much more serious when designing tank- 
age for the low density liquid hydrogen. 

In summary, booster flight loads are seen to be only a part 
of the overall system response, subject to extremely involved 
design compromise and especially dependent upon the ability 
and courage of the control system engineer. With a bold and 
coordinated approach, it appears that problems arising from 
coupling of structural and sloshing modes with the control 
system can be reduced in severity and conquered by sophis- 
ticated design. 


Environmental Conditions 


The forward surfaces of Earth takeoff rockets suffer some 
aerothermal effects which produce relatively modest struc- 
tural design problems. Re-entry problems are in a separate 
category, and have been widely treated, for example, (9). 

The environment of space is a topic which currently is of 
strong interest to structures and materials engineers (10,11). 
The ultra-hard vacuum which applies above a thousand 
miles altitude cannot at present be duplicated on Earth; from 
such experiments as can be conducted it is known that the 
surface properties of materials, which include crack propaga- 
tion and fatigue characteristics, may be strongly affected by 
vacuum. These properties are especially significant in the 
design of pressurized thin shells and may govern structural 
weight and reliability. Ionizing radiation, chiefly that of solar 
origin, will have a detrimental effect on the stability of struc- 
tural and other properties of plastics and organic materials. 
Although incident radiation is itself of negligible consequence 
in the behavior of metals, when in combination with hard 
vacuum it may hasten or influence the changes which occur in 
surface characteristics. Finally, by their combined influence, 
radiation and vacuum may affect surface emissivity, destroy- 
ing radiative thermal control and leading to detrimental 
temperatures. 

Damage by meteoritic particles constitutes a major en- 
vironmental hazard which has received attention with regard 
to identification of the environment (12,13), nature of the 
damage mechanism (14) and protective design (12,15). 
However, none of these studies can yet be applied with engi- 
neering accuracy. The state of the problem is nicely sum- 
marized by Gazley and Bjork (16): “The probability of 
damage to space vehicles....has appreciable uncertainty 
because of the difficulty of estimating meteoroid mass and 
frequency. The results indicate, for example, that a surface 
1 meter square and 1 millimeter thick will be punctured some- 
time between 7 weeks and 620 years.” 

In summary, the interaction of structures and environment 
will exert a major influence on the progress and the cost of 
spaceflight, since it is evident that solution of the environ- 
mental problems could involve a doubling or tripling of the 
dead weight ratio as a result of provision for shielding, meteor 
protection and fail-safe design. As yet, however, the re- 
quirements are undefined, measures to meet them are hardly 
begun, and an effort at quantitative prediction is premature. 


Structural Forms and the Vehicle Design 


Loading Density 


Requirements for the containment of large volumes under 
pressure and the simultaneous transmission of vehicle loads 
establish thin pressurized shells as the dominant structural 
form. Conventional aircraft shell structure of monocoque or 
semimonocoque design represents one method of supporting 
the thrust and bending loads, which is appropriate when the 
internal pressure is a minor factor. In the other extreme, 
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the pressure itself may be used as the medium which supports 
the overall compression and bending, with the shell acting 
chiefly as a membranous container. Successful pressure 
stabilized structures date from the all-metal dirigible of 
1929 (17) and include many contemporary examples. In 
many other cases, internal ullage pressure will support a part 
of the overall axial load or bending moment, leaving the ba'- 
ance—as well as the hoop tension loads—to be carried by the 
shell itself. Which type of structure is most efficient depenc s 
on load, size and shape, all at once. 

In approaching such a design problem, the structures eng)- 
neer often makes use of the “structural index” (18). Th» 
governing structural indexes for the three major load-carryir z 
functions under consideration here are: 


Axial load P/rk? 
Bending 2M /xR* 
Internal pressure p 


The value of the structural index is indicative of the densit y 
of the structure required, and because material properties pe - 
mit only limited freedom in this regard, the structural ind: x 
usually determines the fabricational form which correspon: s 
to minimum weight. The structural indexes selected hei: 
place all three types of loading on a unified basis, namel>, 
equivalent internal pressure. 

The way in which the values of these structural index: s 
vary with the vehicle design parameters of p,, Wer and L/:) 
has been derived in Table 1 for the major types of impose: 
loads. The assumptions on which each is based, also note:| 
in Table 1, follow from considerations of the previous section ; 
these assumptions are simplified but realistic. 


Design Capability 


The price in structural weight which the engineer must pay 
to meet these requirements is expressed graphically by employ- 
ing the analytical theory of optimum design, as in Figs. 1 and 
2. Allowable strength data for conventional stiffened shell 
structure have frequently been cast into such a form [for 
example (19 and 20)]. It is also possible to include data for 
pressure stabilized shell structure as well. The ordinate for 
these curves is a weight parameter p,t/R, where ¢ is the thick- 
ness of a simple unstiffened shell of the same weight (equal to 
the thickness if pressure stabilized), p; is the material density 
and R the shell radius. The abscissae are the structural 
indexes P/rR*? and 2M/zxR*; these parameters will be recog- 
nized as being equal to the stabilization pressure that would 
be required for a shell with no inherent stiffness. 

In constructing these charts the strength data for un- 
pressurized shells were based on (20), which is somewhat 
more conservative than (19); the performance of simple 
monocoque structure and of pressure stabilized shells was 
based on (21 and 22), with the allowable hoop stress limited to 
85 per cent of the material tensile yield stress. The figures 
apply directly to 7075-T6 material, but the effects of other 
materials may be estimated by introducing their properties 
into the ordinates, as is indicated, for example, in (19). The 
curves for pressure stabilized shells of austenitic stainless 
steel of 200,000 psi yield strength are conservatively predicted 
by those for 7075-T6. 

From a cursory study of Figs. 1 and 2 it may be concluded 
that for low values of the structural index (low loads, large 
dimensions) pressure stabilized structure has a significant 
weight advantage. For high values (concentrated loading) 
conventional stiffening is optional. However, when pressure 
is to be carried in addition to axial or bending load, longi- 
tudinal stresses in the conventional stiffened structure are 
relieved, and the effect is to reduce the design value of the 
structural index by the magnitude of the ullage pressure 


(P/R?) ettective = — Pullage 
(2M ettective (2M sctuat — Pullage 
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These effective values of structural index and the correspond- 
ingly lower weight parameters for the stiffened-and-pressur- 
ized design indicate the possibility of a weight advantage over 


a limited range. The resultant design, of course, must not be 
inadequate under the simultaneously applied hoop tension. 
With respect to the gas pressure storage problems, the 
structural indexes contain the storage vessel dimension r 
which can usually be adjusted to exploit material properties 
to the utmost. 


Significance of the Analysis 


‘his presentation of structural performance in terms of 
str:ictural parameters is not meant to serve in place of a stress 
an:lysis. Its purpose is to express the requirements and the 
-aabilities in gross terms, a form more immediately meaning- 
fu! to the designer. 

‘Vhat does the analysis signify with regard to future prob- 
leis of the space vehicle designer? Several conclusions may 
be drawn from Table 1: 

The effect of size on the structural indexes is slow, but a 
de inite increase with size is unavoidable. 

2 The effect of slenderness ratio is mixed, and changes in 
LD will therefore lead to new design compromises. 

> The effects of propellant density and of acceleration 
fa-tor are large, and the reductions which will occur in these 
quantities with the use of liquid hydrogen propellant and 
space-based vehicles will lead to much smaller structural 
indexes. 

Figs. 1 and 2, however, clearly show that, thanks to pressure 
stabilization methods in the usually troublesome low range, 
the designer is well prepared to meet any foreseeable loading 
condition with structures that utilize the available properties 
of materials in an efficient manner. 


Structural Weight and the Vehicle Design 


Bossart (3) has presented a study of the effect of gross size 
on the dead weight ratio of geometrically similar rocket 
vehicles. Without making extensive and detailed studies, 
it is possible to present only a few new considerations, some of 
which were suggested by Bossart, namely, the effects of varia- 
tion in shape or fineness ratio (L/D) of the vehicle and the 
effect of propellant density. It is possible also to resort to 
empirical data, in which case the effect of a continuously 
advancing technology must be recognized. 


Variation of Component Weight With Configuration and 
Size 


Following essentially the same approach as (3) by consider- 
ing structural weight to be comprised of separable provisions 
against several different major types of loading, and making 
simplified statements concerning the effect of basic vehicle 
proportions and sizes on these structural requirements, the 
structural weight ratios of Table 1 may be derived. The 
propellant mass is seen to be the natural normalizing factor 
for most of the weight allowances (a few depend on gross 
mass, but the propellant loading factor is relatively inflexible 
for a given mission). In general, each weight ratio is roughly 
proportional to the material strength-density ratio, the load- 
ing severity and some function of the vehicle parameters W p, 
L/D and Pp- 

Assuming similitude of loads, materials, design efficiency 
and fabrication, the effects of vehicle design are seen to be as 
follows: 

1 The weight factors allocated for all structure designed 
hy longitudinal thrust, and that designed by bending under 
maximum control system deflection, vary as the one-third 
power of the stage propellant mass Wp. None of the other 
structural weights is influenced directly by size increase. 

2 A decrease in the fineness ratio causes large reductions 
in the weight allocation for material designed by bending, and 
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Fig. 1 Structural efficiency of shells in bending: Optimum 
pressure stabilized structure compared with optimum conven- 


tional aircraft construction 
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Fig. 2 Structural efficiency of shells in compression: Optimum 


pressure stabilized structure compared with optimum conven- 
tional aircraft construction 


somewhat smaller reductions in that applying to longitudinal 
thrust structure, tank bottoms excepted. Tank ends suffer 
with a decrease in fineness ratio. A design compromise 
appears possible which would permit offsetting the weight 
penalties of increased gross size with a reduction in fineness 
ratio. 

3 A decrease in propellant density affects all weight factors 
adversely, except pressure storage vessels. Structures critical 
under tank ullage pressure and aerodynamic loading are 
affected inversely as the first power of the propellant density, 
whereas factors for thrust and hydrostatic structure vary 
only inversely as the one-third power; therefore liquid hydro- 
gen tankage should be designed for low internal pressure and 
small fineness ratio. 

4 A decrease in acceleration factor will bring numerous 
decreases in component weights. 
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Empirical Variation in Body Weight With Size 


‘The present open-literature picture of the state of the art in 
structural design for liquid propellant rockets is summarized 
by Fig. 3, where the body weight ratio is plotted against stage 
propellant weight Wp for a number of contemporary vehicles 
(23). [The body weight ratio is defined as the ratio of 
structure-plus-engine mass to propellant mass (W, + W~,)/ 
W », and does not include residual propellants or fixed-weight 
items, such as guidance and control equipment.] It is evident 
that the larger size vehicles offer a great advantage in struc- 
tural efficiency, and that the trend appears to continue up to 
the largest devices currently visualized. 

Chis observed trend is in contrast with (3 and 24), which 
in licate that, as size increases above present-day applica- 
tions, the structures and engine weight fractions will increase. 
Conservative predictions must assume uniformity of design 
philosophy and geometric similarity overall or in detail, 
whereas in actuality the structure will vary in nature and in 
shape to suit an optimum overall compromise with all aspects 
of the design problem. Structures engineers have never re- 
signed themselves to the law of squares and cubes. More- 
over, such conservative predictions ignore a very real influ- 
ence which the more generous research and development 
budget for a large device has upon the mean level of efficiency 
ol tained in the end product. 

An additional quantity affecting the nature of this chart is 
that, in general, the larger devices have been developed more 
recently in history and benefit from advances in the general 
stute of the art. It is dangerous to extrapolate technological 
advance but unrealistic to ignore it. 


Effect of Structural Weight on Vehicle Design 


The importance of large size booster rockets in the space 
effort is now well recognized. Fig. 3 indicates that ambitious 
projects will also be rewarded with exceptional structural 
efficiency. This effect would be extremely significant to the 
overall performance of a device such as the V-2, which had a 
payload ratio of 0.08 and a body weight ratio of 0.30. But in 
contemporary booster designs the values may be reversed, 
and then the improvement in payload owing to structural 
weight saving, particularly weight saving in the booster, is 
marginal. Instead of exploiting structural efficiency directly, 
it may be more attractive to effect a tradeoff with such other 
desirables as reliability, economy of construction and recovery 
provisions. 

One such alternative may be a reduction in the number of 
stages to perform a given mission. The primary function of 
staging, with its almost direct increase in complexity, cost and 
unreliability, is to avoid the penalty which is incurred by only 
a few per cent of structural weight. A representative ex- 
ample is summarized in Table 2, where for a typical mission 
the ratio of gross takeoff mass to payload mass is tabulated for 
three-stage and two-stage satellite launching rockets of several 
different sizes, assuming body weight ratios given in Fig. 3, 
with appropriate allowance for residual and fixed mass as 
indicated. The results, although idealistic, demonstrate that 
when gross size is increased by an order of magnitude, a 
mission formerly requiring three stages can probably be done 
with only two. 

Low structural mass fractions, furthermore, result in a 
notable reduction in sensitivity of the overall performance to 
staging proportions, and therefore add versatility to the 
attributes of very large boosters. Fig. 4 illustrates this effect 
in a two-stage rocket booster, designed to furnish an ideal 
velocity of 20,000 fps—suitable for launching a third-stage 
satellite or escape vehicle. For several different dead weight 
ratios, overall payload ratio is plotted against the staging 
proportions, the latter being measured by the first-stage mass 
ratio MR, normalized with respect to the mean mass ratio 


VMR, = V Mh MR. (A value of 1.0 for the parameter 
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Fig. 3 Effect of vehicle size on body weight for Earth takeoff 
rockets using conventional liquid propellants 
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Fig.4 Effect of dead weight ratio on staging optimization for two 
lower stages of a satellite launcher 


Table 2 Two-stage vs. three-stage satellite launcher 
performance 


(Based on dead weight estimates 6 vs. W > of Fig. 3) 
Initial mass at takeoff W;, 108 108 107 


Initial mass/payload in orbit, 


three-stage vehicle 90 40 35 
two-stage vehicle 140 50 40 
Assumptions 


= 300 sec 
I2 = I,3 = 320 sec 7 
(vi)o = (va + v2) = 32,000 fps 
where for each stage 
A}=1- (1 -—1/MR)(1 + 4) 
vi = Igo ln MR 


Overall payload ratio X» = A;AsA3; Maximized by adjusting 
staging proportions, 


MR; / MR, means that the mass ratios of first and second 
stages are identical.) Propulsion performance used in the 
examples is roughly representative of combinations of Lox/ 
RP-1 and Lox/hydrogen systems, with allowance for altitude 
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Fig. 5 Family of large booster structures (designed for thrust = 2W>p approximately 10 deg misaligned) 

Clearly, first-stage booster vehicles of low structural weight Representative Booster Designs Pa 

may be utilized to carry a wide range of off-optimum payloads ; 

with small sacrifice in performance, whereas a large dead As an example of the effects of size on structure and 0 
weight fraction will amount to a restriction in utility as well structural problems on configuration, a series of cursory de- 
as a direct performance penalty. signs for single-shell booster vehicles was made, covering th: 


range from 10°- to 10®-lb propellant mass. Fig. 5 presents the 
results of these design studies and indicates trends which maj 


System Design vs. Structures 
“i . be expected for such vehicles. Each design is the result of a 


The evidence is strong that the dead weight fraction for number of trials and represents a degree of refinement, insofa1 
liquid propellant rockets should go down, not up, with size as the depth and limitations of the study permit. 
increase. Dead weight is, currently, the primary reason for Conservative methods and assumptions were followed in 
staging and for carefully optimized staging proportions. Ina developing these designs, and the stipulated loadings were 
space systems design the structures contribution must be met with factors of safety of 1.25 and 1.50, respectively, for 
recognized and exploited, and this means much larger vehicle unpressurized and pressure-carrying structures. This factor 
sizes and fewer units than are usually considered feasible for of safety requirement, together with the conservatism of 
chemical rockets. Figs. 1 and 2 regarding the weight of pressure stabilized shells, 
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Table 3 Structural weight estimates for booster designs of Figs. 5 and 6 


Factor of safety 1.50 on pressure loads; 1.25 on other loads. 
bending stress corresponds to about 10-deg thrust misalignment. 


Propellant mass Wp, lbm 105 10— — 10’ 108 
Propellant Lox/RP-1 Lox/RP-1——————._ Lox/hydrogen Lox/RP-1 Lox/RP-1 
Booster length, ft 57 94 75 65 109 125 225 
Booster diameter, ft 83 18 22 26 30 55-35 135-80 
Ullage pressure, psi 45-45 20-20 20-20 45-45 34-34 0-60 0-60 
Applied axial load for shell, lb 2 X 105 2 X 106 2X 10° 2 X 108 2 X 10° 2 X 10° 2 X 108 
Applied bending stress, % of A 
applied compression stress 85 75 ‘ 75 —6% 75 50 33 
Design axial load equivalent to 
combined axial and bending 7 5 
stress 3.28 3.18 10° 3.43 X10 3.18 K 10 10" 275K 32:5. xX 
Tankage weight, Ibm 
Longitudinal shells, axial 
plus pressure loads 670 15,400 10,500 3800 6900 107,000 1,550,000 
Tank bottoms, pressure loads 
only 296 3370 4990 8200 11,600 81,000 1,234,000 
Domes, ullage pressure only a 225 760 1130 1820 1760 Mae sist 
Other structure weight, Ibm 
Longitudinal shells: Inter- _ 
stage, aft struct. S10 5380 6400 7750 8100 27,700 87,000 
Engine mount and skirt 340 ~=~—«~*:70 2320 2820 3050 15,000 80,000 
Miscellaneous plus 5% con- t ~ 
tingency 300 2450 2400 2490 3250 21,200 240,000 
Ullage pressure storage system —_ 
weight, Ibm 900  —-2500 2500 5500 13,600 19,800 350,000 
Total structural weight, Ibm 3541 31,630 30,240 32,380 48 , 260 271,700 3,541,000 


Longitudinal loads based on thrust F = 2Wp; superimposed 


predisposes the larger designs to stiffened-shell structure. 

Weight calculations based on these studies support Fig. 3 
out to 10’ lbm. For the 108-lbm design the structural weight 
fraction is greater than that for the 10’-lbm vehicle. The 
difference is not so large as to discourage more thorough 
study and consideration. 

A typical crossplot is shown in Fig. 6, where a number of 
10®-Ibm boosters of different L/D are compared. The 
weight breakdown for these designs, given in Table 3, sub- 
stantiates the detail trends predicted by Table 2. The over- 
all result, however, indicates a surprisingly slow variation in 
total weight fraction with L/D for the range investigated. 
The other studies did not permit a similar observation, the 
10°-lbm booster design being constrained by minimum gage 
considerations, and the larger vehicles by engine size and 
tankage shape problems. 

Fig. 6 also presents, for comparison, a design using liquid 
hydrogen fuel, at a compromise oxidizer to fuel mixture ratio 
of 5. About three times larger in volume than the oxygen- 
kerosene rocket of the same mass, the structure shows 60 
per cent increase in weight. This result is in accordance with 
Table 1, considering that most of the structure carries longi- 
tudinal and hydrostatic loads and its weight would vary as 
Vv, Pp if the material could be utilized with equal efficiency. 
Pressure stabilization makes this roughly possible. 

In a study of this type it is not possible to anticipate the 
ingenuity of the design engineer, who by combining functions 
and rearranging configurations may gain a notable advantage. 
It is presumed that many technological advances will be 
achieved, however, for the larger devices in the family: 
Thrust damping of propellant sloshing, ullage pressure con- 
trol, pressurized interstage structure, double-wall tank bulk- 
heads with internal pressure, hot gas pressurization systems, 
and the plug nozzle, to transmit engine thrust directly into 
the tankage, eliminating the tortuous load path that is toler- 


ted in contemporary vehicles. = 
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Conclusions 


It appears that, at present, structures place no bounds on 
the size of rocket vehicles, and in fact favor much larger 
devices than are presently contemplated in our space pro- 
gram. Dynamic problems arising from the interaction of 
structural deformation, fuel sloshing modes and the control 
system will probably constitute the area of greatest technical 
difficulty associated with new, unusual and _ structurally 
efficient designs. From the crude analysis presented here, it 
appears that on the whole these problems may be relieved 
with increase in size. Despite this possibility, an advance is 
needed in controls systems to provide isolation and active 
damping measures. 

Many advances have been made since the day of von 
Pirquet. Mammoth booster vehicles, such as those studied 
here could establish large lunar and interplanetary missions 
not only without requiring orbital assembly, but probably at 
reduced total cost, and certainly much earlier in historical 
time. 

The structural problems of large space vehicles employing 
thin wall tankage will be much the same as those for today’s 
boosters, differing chiefly in the values of the design param- 
eters. Efficient methods of approaching these new design 
problems have already been established, although the degree 
of efficiency can undoubtedly be improved. There remain a 
number of questions concerning environmental effects, how- 
ever, which cannot be treated quantitatively at this time. 
The recent NASA Technical Note D-518 (10) describes in 
detail what would be necessary to approach these problems 
with assurance and competence. 

Actually it is not necessary that the structures engineer be 
well equipped for this task. The alternative is greater struc- 
tural weight in order to obtain a prescribed confidence level. 
For a given mission, this requires either more complex staging, 
a design approach which must face up to more severe relia- 
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bility limitations and mission cost, or improved propulsion 
performance, which again can be measured in cost and time. 
The possibility of a bold structures approach contributing to 
the overall performance of new design is in danger of being 
overlooked under the pressure of meeting contemporary 
problems with conventional methods. In brief, structures 
technology is not a stationary science, and the conquest of 
space can progress better, faster and more economically when 
founded on coordinated advances in all fields. 
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Nomenclature 

Cma = slope of aerodynamic moment coefficient 

D diameter of vehicle 

E = modulus of elasticity oy 

F, = allowable stress ae 

I, = specific impulse, lbf-sec/lbm he 
Iven = mass moment inertia vehicle, slug-ft? in 


Tsnei) = section moment inertia shell, 


K ~~ = bulk modulus of fluid, psf 

= length of vehicle 
= bending moment 

MR = mass ratio = W:i/W, 
P = longitudinal thrust load, lbf 
R = radius of tank, in. 
7, = natural period, sec 
W; = initial gross vehicle mass = W, + We, Ibm ih mi 
We = useful propellant mass, hen 

W, = burnout mass = W, + Wa, lbm 7 =n 
= payload mass, lbm 

Wa = “dead weight” mass = W, + We + W,, lbm 

W, = mass of vehicle structure including tankage, lbm 

We = engine mass, lbm 

W, = residual and fixed mass, lbm 

f = frequency, eps 

g = gravitational conversion constant 

h = height of liquid level, ft 

k = constant of proportionality 

n = acceleration factor 

Pp = internal pressure, psi 

q = dynamic pressure = 
r = radius of pressure sphere : 

t = skin thickness, in. 

vi = ideal velocity = (J.g) In MR, fps 

f.s. = factor of safety 


= angle of attack 
= swivel engine control 


Rocket,” 


Pp 
ps = density of structural material, lbm/in.* 
o = fiber stress 
T = time, sec a 
6 = dead weight fraction, Wa/Wp 
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tended to more general situations. 


3 THE following, some of the basic phenomena of satellite 
:otion are discussed in the most direct and unsophisti- 
catcd way that seems possible. The derivations are really 
the outcome of a quest for independent checks of the results 
obtainable from the well-established and refined methods of 
astronomical perturbation analysis. Naturally, an approach 
from first principles, i.e., energy and momentum conditions, 
seemed desirable. It is felt that these developments may be 
of tutorial value to the non-expert who desires a basic under- 
standing of satellite orbit perturbations (on the other hand 
this “return to nature” may have some amusement value for 
the specialist in dynamical astronomy). The physical 
phenomena treated here are: The distortion of an orbit 
within its own plane resulting from extraterrestrial gravita- 
tion, apsidal precession and buildup of eccentricity resulting 
from radiation pressure, apsidal precession and changes in 
epoch resulting from the equatorial bulge, and nodal preces- 
sions resulting from oblateness, extraterrestrial gravitation 
and radiation pressure. The geometric relations between 
the orbit and the source of perturbation have been specialized 
for maximum simplicity in each case. The quantitative re- 
sults are precisely correct in all but one example, where a very 
good approximation of the true rate of precession is obtained. 


In-the-Plane Perturbations of the Orbit 


For the demonstrations of in-the-plane perturbations of the 
satellite orbit, we start with the simplest shape—a circular 
orbit. A convenient coordinate system by which the satel- 
lite’s departure from the unperturbed orbit may be described 
are the £,n axes of Fig. 1. Their origin lies at a point such as 
0’ which would be the vehicle’s instantaneous position in the 
absence of perturbations. The £ system then is moving 
with 0’ while the z,y system, with origin at Z, may be taken 
as an inertial one. The instantaneous radial and tangential 
velocity components of the satellite shall be denoted by U, 
and U;; the perturbative forces in the £ and 2,y system 
shall be the symbol F with suitable subscripts. (See Fig. 1.) 
We may then derive the equations of motion of the vehicle 
in the £,» system from the fact that the time derivative of its 
energy K equals the rate at which work is done on it by the 
perturbative forces, and the time derivative of its angular 
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* For treatments of this kind, the reader is referred to such 
well-known texts as references (1 and 2). Since the present dis- 
cussion is not an attempt to emulate the classical ip our 
somewhat unconventional notation 
some astronomers. 
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Fundamental Satellite Perturbations 


: A brief exposition, from an elementary point of view, is given for some basic satellite perturba- 
The treatment serves to generate independent checks of formal results from more elaborate 
analyses and to establish a direct relation of these results to physical principles. 
to in-plane and precessional perturbations of circular orbits resulting from oblateness, luni-solar 
gravitation and radiation pressure. T 
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he derivations have quantitative validity and may be ex- 
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momentum A _ E equals the moment of the perturbative 


ee 


forces about #. Thus 
(1] 
dA/dt = Fire +’) —Fm [2] 


where 7, is the radius of the circular orbit. | Clearly 


GEm 

+ + n')(re + (¢’ n)n] 


where the primes indicate differentiation with respect to 8, 
dots denote time derivatives, and GE is the universal gravita- 
tional constant times Earth’s mass. Substituting Equations 


[4] 


[3 and 4] into Equations [1 and 2], we find after some simpli- 


v 


fication 


GEm 


+ — + + 9”) 
nt” + én”) = Fire + E+ + — 


and 
GEm 


r3 


[(re + + 0”) + (20 — = 
Fi(r, + §) — Fin 


[6] 


If we retain terms of O(r,£) in the brackets of the left-hand 
sides and O(r,) in the parentheses of the right-hand sides, both 
equations yield 

+ 9” = (r3/GEm)F;, [7] 
Going to 0(£?) on the left and 0(£) on the right, we may elim- 
inate between Equations [5 and 6] to obtain 


— — 3§ = GEm)F, [8] 


The solution for the homogeneous form of Equations [7 and 
8] may be easily found as 


(r,3/ 


is +38] + 2 [9] 
2, 1 


2? te + 33 | sin 6 + cos @ [10] 
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This may be used in a straightforward manner to generate par- 
ticular solutions for nonvanishing right-hand sides, which are 
really the cases that interest us here. 

As mentioned before, our first example will be the perturba- 
tion resulting from extraterrestrial gravitation. We locate a 
disturbing body with mass P on the positive y axis (Fig. 1) so 
that 
F,= sin 0 

TEP TEP 

where rzp is the distance from E to P. The forces F, and F, 
are the components of the vector difference of perturbative at- 
traction by P on the satellite and on Earth. Transforming 
this to F; and F,, and substituting these expressions into 
Equations [7 and 8], we find the solutions 


F, = — 


= = (1/4)g [80 + 11 sin 26] 


where g = Pr,4/2Ergp*. The first term in each of these ex- 
pressions is secular, whereas the second one is periodic in 
nature. The periodic term in ¢ shows that the orbit develops 
a major axis in the z direction and a minor axis in the y direc- 
tion. The periodic term in 7 expresses the angular anomalies 
which are caused by the varying distance from E to the 
satellite, r, + & in order to fulfill the conditions on A and 
K. This treatment may be easily extended to allow for a 
moving body P, which, however, adds little to the demonstra- 
tion of fundamental principles. If we apply these results to 
lunar perturbations on a 24-hr satellite, we find that the 
maximum amplitudes of this effect are approximately 3000 ft 
in 7 and 2200 ft in ¢ Both displacements are negligible 
from a practical point of view, whereas the effects of inaccu- 
rate launch conditions may be much more significant, as cal- 
culated from Equations [9 and 10]. 

In order to study the effects of radiation pressure, we con- 
sider the same configuration as in Fig. 1. The perturbative 
forces are then 


F, = —fm/r*gp F,=0 


where £ is a coefficient which accounts for radiation intensity 
and shape, size, surface to mass ratio, and reflectivity of the 
satellite. From this we note that 


which when substituted into Equations [7 and 8] yield 
= g [63 cos 6 — (7/2) sin 6] [13] 
n = — g(66sin 6 + 11 cos 6] [14] 


where g = 6r,3/2GEr pp’. 

The most significant features about this result are the 
periodic terms with increasing amplitudes in both displace- 
ments. They are not a surprise if one notes that the excita- 
tion given by F, and F; isin resonance with the free oscillations 
of the satellite about the point 0’ as shown by Equations [9 
and 10]. Of course one then expects a response whose am- 
plitude grows without limit. It may be shown that this fea- 
ture is indeed not far from physical reality if the effects of 
other perturbations are not considered simultaneously with 
radiation pressure. For the time being, one can deduce the 
following useful results from this analysis. = 

The condition d¢/dé@ = 0 reads 


6=—(1/6)cot@ (15) 


where the pair of roots for 6 in each interval of 2x gives the 
locations of perigee and apogee in that particular pass around 
the orbit. Inspection of Equation [15] shows that the two 


3 The satellite is assumed to be an isotropic reflector. i 
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= 


for large numbers of revolutions. 
of roots for w, reads 


+3134 


where the terms n(__)/r, are small corrections and cf 


lim w, = (2n — 
n is the number of the current revolution. A corresponding 
sequence of roots may be written for we, the apogee locaticn, 
with 
lim = 
@ 

A physical interpretation of the way in which the apsides «re 
formed by radiation pressure may be given in the following 
manner. The satellite is being accelerated tangentially »y 
F, during the approximate interval 7/2 <6 < x. In orcer 
to balance the additional centrifugal force resulting from 
the increased circumferential speed, a decrease in r is 
required to enhance the gravitational attraction. Thus 
perigee is formed in the neighborhood of @ = x. The oppo- 
site mechanism acts to generate apogee near 0 = 27. 

Now let us define the average value of the semimajor axis 
for one revolution as 


@=rt (1/2) (Emax + &min) 


and let us assume for convenience that wg = w, + 7 in each 
revolution, then 


a ~r, — 91(3/2) cos wy [16] 


which proves that the orbit size is practically unaffected by 


radiation pressure. For the eccentricity, averaged over 27, 
we get however 


1 g 
é= (max — ~ COS wy + 7 sin 
[17] 
This indicates a secular buildup of eccentricity at the rate 
é’ (—3g /r,) cos We 


Again we should emphasize that the validity of this result 
needs to be examined in the light of the various factors that 
enter a complete solution of the problem. They are the vari- 
ation of the gravitational potential gradient for large depar- 
tures from the circular orbit, variable radiation intensity as a 
function of @ resulting from Earth’s shadow and albedo,‘ and 
the interaction of perturbations from physical sources other 
than radiation with the effects of the latter. Clearly sucha 
complete treatment is beyond the scope of the basic demon- 
strations intended here and is given elsewhere. [See refer- 
ence (3).] 

Our last example of in-the-plane perturbations of near- 
circular orbits concerns the precession of apsides for a slightly 
elliptic orbit owing to Earth’s oblateness. If the oblateness 
term were absent from Earth’s potential, then a slightly 
elliptic orbit may be described by Equation [10] as resulting 
from slightly noncircular launch conditions. Here the circu- 
lar frequency of the periodic terms for £ and » is that of the 
orbital motion, and the axes of the eccentric orbit would con- 
sequently remain stationary in the z,y plane. If, however, 
the restoring forces in the satellite’s oscillations about the 
circular orbit are slightly altered by the oblateness term, 
then a corresponding frequency change occurs in £ and », and 
the latter are no longer in tune with the orbital motion of 0’. 
Consequently, a precession of the apsides results. Our first 


‘ This feature turns out to have —_ a minor effect on the 
results [15 to 17]. 
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«Fig. 1 In-plane perturbation 


task then consists of modifying Equations [7 and 8] to account 
for : nonspherical Earth. This we may do by using the ap- 
pro:ch which led to Equations [5 and 6] but replacing the 
last term in Equation [3] by a more accurate representation 
of the geopotential. This potential function is usually as- 
sumed to be rotationally symmetric about Earth’s polar 
axis, here the z-axis, and thus lends itself to expansion in 
terms of Legendre polynomials in sin ¢ where ¢ is the geocen- 
tric latitude. One then finds 


where 


2/r = sing 

Te = 3963.247 miles, Earth’s equatorial radius 
J = 1.638 10-3, the coefficient for P» 

D = 1.070 10-5, the coefficient for Ps 


For the present discussion we may safely drop D from Equa- 
tion [18], so that for z = 0 the radial equilibrium condition 
of the nominal orbit can be written as 

Vie = GEm (1 + 


c 


Hence 


which indicates the slight correction to be made in 6, the 
angular rate of 0’, as compared to motion around a spherical 
Earth. Writing the time derivative of Equation [18] and 
using r?> = (re + £)? + 7? to write binomial expansions for 
and 1/r5, one obtains 


to 0(é*) in the brackets. Using Equation [19] this becomes 


= mé? [ret — 2(1 + + [20] 


where X = Jrg?/(r2 + Jrg*). Substituting Equation [20] 
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as the potential term in the time derivative of Equation [3], 
we arrive at the following modified version of Equation [8] 


for the homogeneous case 


£” — 2n’ — (8 + 24) = 0 


= 


which of course reduces to [8] if J = 0. The homogeneous 
form of Equation [7] is unaffected by the change in potential. 
The solution of this pair of equations is a 


= b, sin + be cos x6 + bs 


2b 
n= by — 5 (3 + — sin + cos 
kK K 


where «x = 1 — A to O(A), and b;.2.3.4 are integration constants. 

Since the circular frequency of and 7 is less than the orbital 
frequency, a forward motion of the apsides occurs which cor- 
responds to the average rate 


= = bJrg?/(r2 + Ire?) [22] 


This result differs from the first-order approximation by 
more elaborate techniques only by the presence of Jr,” in the 
denominator. 

For a 6000-mile orbit and the constants that were given with 
Equation [18], we obtain w = 0.479 10~ rad per sec, and thus 
an apsidal period of 152 days. Finally, the decrement of 
epoch during a specified revolution amounts to 


[n, (0) — »,_,(0,)] [23] 


where 6; marks the point whose time of passage is used as the 
epoch, and n identifies the orbital pass that is being com- 
pleted. 


Perturbations of the Orbit Plane ae 5 


We will now discuss the physical effects which force a 
satellite to depart from its original orbital plane, or in other 
words, which change the orientation of this plane. We now 
tilt the circular orbit about the z-axis and denote the angle of 
inclination to the x,y plane by 7 (see Fig. 2). The ascending 
node is N. Applying the right-hand rule to the rotation of 
the satellite in its orbit, the angular momentum vector A of 
this motion forms the angle 7 with the positive z-axis. We 
now return to Equation [2] in its vectorial form 


dA/dt = M 
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Fig. 2 Perturbation of the orbit plane 
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where M is a moment that may be exerted on the orbit by 


perturbative forces. If these perturbations result from the 
equatorial bulge, we observe from symmetry that the average 
moment over one revolution lies along the z-axis. It shall 
be denoted by ,. Ifa positive M; is taken to point in the 
positive x direction, we note by the preceding equation that 
it produces a rotation of the momentum vector about the z- 
axis at a rate Q, = M,/A sin i while the magnitudes of the 
momentum vector and the inclination remain unchanged. 
The positive sense of Q, is defined again by the right-hand rule 
about z. The nodal points of the orbit will then precess in 
the x,y plane at the rate Q = Q,. 

Since 6 is constant around the orbit, a time average of 
| M | is equivalent to an average against 6. As a result of 
Earth’s oblateness we have 


F, = —V, = —(2GEmJr,?/r,') sin @ sin 7 
. GEmJrz? . .. 
i, = 2i sin? 6d@ = — sin 21 
2ar,3 0 2r¢ 


(The contribution of F, to this moment can be shown to 
vanish.) Recalling that 

A = 6r2m = (GEr,)'?m 
we find 


Since Q is negative, we speak of a regression of the nodes. 
This result agrees exactly with the expression which can be 
obtained by more advanced methods. If r, = 5000 miles 
and i = 51 deg, we find Q to be equal to —5.80 107 rad per 
sec or 0.24 deg per revolution (where the anomalistic period 
is 7180 sec). The nodal period amounts to 125 days. 

Finally, we examine the nodal precessions resulting from 
gravitation and radiation from a distant body P which we 
locate again on the positive y-axis. A perturbative moment 
M, will be exerted on the orbit, resulting in a precession 


Q, = M,/A cosi [25] 


which in this case takes place about the y-axis because the 
moment vector, initially M,, remains always in the z,z plane 
as the orbit is perturbed. Inspection of the geometry at the 
node N shows that for small precessional angles we have 


and 


= Q,coti = sini [26] 


For the case of cipnacesbe interference, we note that 


7 di/dt = 0 


F, makes no contribution to M,. Then 


M, (Fyr, sin 6 sin t — F,r, sin 6 cos 
0 


3 sin 2ir,2GPm 
4r EP® 


With Equation [26] this results in 
— (3P/2Er gp*)(GEr,3)'? cos 


ll 


For a 24-hr satellite with i = 20 deg we thus find Q = 
10~° rad per see owing to lunar gravitation. 


- magnus force on spinning satellites. 


Turning finally to the precessional effects of radiation pres- 
sure, we note that 7, = 0 for a uniformly illuminated orbit. 
Since this case is not at all realistic, we add the factor (1 + 
sin @) to the radiation pressure, indiating that the peak intcn- 
sity occurs on the side of the orbit toward the illuminating 
body. Thus 


TEP” 


The sine term in F, will now generate a perturbing mom: nt 
on the orbit which is 


M, = (8mr,/2rgp”) sin t 


Using Equation [26] : | 


1/2 
[0] 
2r GE 
It is interesting to note that this precession is independ: nt 
of 7. For a 100-ft mylar balloon with a metallized skin of 
0.0005-in. thickness, the nodal precession of a 24-hr or it 
would amount to 0.6325 10~7 rad per sec. A more compl: te 
discussion of the cumulative precession due to radiation prvs- 
sure is given in reference 3. 


The elementary techniques discussed may be readily applied 
to other physical effects, such as atmospheric drag and the 
It is hardly necessary 
to point out that the simple approaches which were given 
here could, with due perseverance, be extended to the more 
complicated geometry of more general situations, although 
they do not compete in elegance and power with the estab- 


lished methods of dynamical astronomy. They have, how- 
ever, served us well in providing spot checks by which most 
perturbation phenomena may be traced back to elementary 
principles. 
Nomenclature 
m = satellite mass ae 
GE = mass of Earth times universal gravitational constant 
GP = mass of planet times universal gravitational constant 
x,y,z = geocentric Cartesian coordinates 
€.7 = moving coordinates relative to nominal satellite posi- 
tion 
= distance Earth to satellite 
TEP = distance Earth to planet 
Te = radius of nominal circular orbit 
0 = nominal angular position of satellite 
K = total energy of satellite 
A = angular momentum of satellite | 
= perturbing forces acting on satellite 
= radial and circumferential velocity of 
satellite 
TE = Earth’s equatorial radius | 
V = Earth’s gravitational potential | 
J,D Earth’s oblateness parameters 
#g.@ = locations of perigee and apogee 
@ = apsidal rate of precession 
Q = nodal rate of precession 
References 7 
1 Moulton, F. R., ‘Introduction to Celestial Mechanics,” MaeMillan 
Co., N. Y., 1958. 
2 Smart, W. M., ‘‘Celestial Mechanics,’’ Longmans, Green and (0., I 
N.Y., 1953. 
3 Geyling, F. T., ‘Satellite Perturbations from Extra-Terrestrial Gravi- t 
— and Radiation Pressure,”’ J. Franklin Inst., vol. 269, May 1960. p. . 
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Qotimum Orbital Transfer by Impulses 


Polytechnic Institute of Brooklyn 
Brooklyn, N. Y. 


The minimum impulse transfer from a given initial orbit to the final orbit of given geometry ac- 
complished by multiple or single impulse is considered. It is shown that the total impulse will be 


am a minimum if the two orbits and the transfer trajectory are co-planar. If there is only a single im- 

[. pulse available, the orbital transfer can be accomplished provided that the final orbit can intersect 

Lae 7 the initial orbit at some orientations. The procedure of finding the optimum orientation is pre- 

e — sented. With two impulses available, the orbital transfer can be accomplished without any limi- 
- oe : tation on the geometries of the initial and the final orbit. It is shown that the sum of the two im- 


pulses will be a minimum if the major axes of the two orbits and the transfer ellipse coincide with 
each other and are oriented in the same sense, and the apogee of the transfer ellipse coincides 
with that of the orbit with longer apogee distance while the perigee of the transfer orbit 


coincides with that of the other orbit. a nn -— 


HE PROBLEM of transferring a vehicle from one orbit Proof: 
into another about the same attracting body under the 


For a given system, or a given value of the sta- 
tionary mass M at the focus, the areal velocity h is constant 
for each orbit. 


optimum condition has been investigated by Hohmann (1),? 
Lawden (2,3), Munick (4) and others [see references in (5) ]. 
In their investigations, the orbits are co-planar and explicit, 
or numerical results have been obtained only for the special 
case of two circular orbits (1) and the case of two elliptical 
orbits with special geometry and/or orientation (2-4). 

In the problem of re-entry of a vehicle from an initial orbit 
to Earth, it is desirable to transfer the vehicle into a final 
orbit which will enter the “limit of atmosphere” near its 
perigee, so that atmospheric drag can be employed to reduce 
the velocity of the vehicle. The geometries of an orbit are 
defined by its eccentricity e and its semilatus rectum p. 
For the initial orbit, they are specified. For the final orbit, 
they are determined by the entry velocity and angle at the 
limit of atmosphere. The orientation of the final orbit with 
respect to the initial orbit is not specified. It should be such 
that the orbital transfer takes place with minimum fuel con- 
sumption. 

In the present paper, it is assumed that the initial orbit 
and the transfer trajectory lie outside the limit of atmosphere 
and that the orbital transfer is accomplished by several im- 
pulses. The transfer trajectory will be composed of several 
elliptical ares which together with the initial and the final orbit 
possess a common focus. The change in velocity at the 
intersection of two elliptical ares is proportional to the im- 
pulse which represents an integral of the thrust of a rocket 
engine during a very short firing period. If the rocket en- 
gines are similar, the requirement of minimum fuel consump- 
tion is equivalent to that of minimum impulses required for 
the transfer, which in turn is equivalent to that minimum 
change in velocity. 

The investigation begins with the proof that the optimum 
orbital transfer can be reduced to a planar problem. 


Justification of Reduction to Planar Problem 


To prove that for the optimum transfer the two orbits 
and the transfer trajectory should be co-planar, it is shown 
first that for any orbital transfer at the intersection of two 
non co-planar orbits there exists a corresponding co-planar 
transfer with less impulse required. : 


Received Jan. 20, 1960. 
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The equation of the orbit and the radial 
and circumferential components, v and u, of the orbiting 
velocity are (6) 


[1] 
+ e cos [2a] 
v=r=-—h (*) = [2b] 
di \r Pp 
h? h? 
WG [3] 
where 
rf = distance from the common focus ' 
= true anomaly > 


G = universal gravitational constant 


u and v depend only on the angular distance from the peri- 
helion. 

As shown in Fig. 1, orbits 1 and 2 intersect each other at 
point K. Therefore SK is the line of intersection of the plane 
of orbit 1 and that of orbit 2. The angle between these two 
planes is denoted by ¥y. The radial distance of AK from the 
focus S is r* with 

1 1 1 

= = — (1 + cos = — (1 + cos [4] 

Pr Pe 
where the starred quantities refer to those at the point A, 
and subscripts 1 and 2 indicate quantities associated with 
orbits 1 and 2, respectively. For both orbits, the directions 
of their radial velocity components at the point K are the 
same, whereas the directions of their circumferential velocity 
components lie in the plane normal to the radial direction SK. 

Since the direction of the circumferential velocity lies in 
the plane of the corresponding orbit and normal to the line of 
intersection of SK, the angle between the velocities uw. and wm 
is y. The change in velocity for the orbital transfer at that 
orientation is 


V — vo*)? + — cos + us*? [5] 
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If the plane of orbit 2 is rotated about the line SK while the 
plane of orbit 1 is stationary, point K remains at the point 
of intersection of the two orbits. The focus of both orbits 
remains at S; and the value of 3,* and 3,*, and hence those 
of u,*, vi*, w* and v2*, are unchanged. In the process of 
rotation, only y is changed. If w*, w,* have the same sign, 
i.e., if the two orbits are in the same sense, AV attains its 
minimum value at y = 0;i.e., the two orbits are co-planar. 

If u,*, ue* have opposite sign, AV attains its minimum value 
at y = x. Again the two orbits are co-planar. This con- 
cludes the proof. 

If the transfer from the initial to the final orbit is accom- 
plished by several impulses and the transfer trajectory is 
composed of several elliptical arcs, then, according to the 


Fig. 2 Orbits 1 and 2 and the transfer orbit 3 


Wigs 


the optimum transfer takes place when |AV]| or (AJ)? 
_isa minimum, i.e. 


and - 


With these two relationships, it can be shown that 


i For the upper half of orbit 1,0 < 3, < x; therefore the first 

"derivative of AV is always positive, and the optimum transfer 
takes place at an angle 8,* less than 3;. For the exceptional 


preceding proof, the two elliptical arcs associated with the 
same impulse should be co-planar. Therefore, the following 
conclusion can be drawn: For the optimum orbital transfer, 
the initial orbit, the final orbit and the transfer trajectory should 
be co-planar. 

The next step is to find the optimum transfer trajectory 
and the relative orientation of two co-planar orbits. 


Optimum Transfer by Single Impulse 


With single impulse the transfer trajectory reduces to t/ie 
point of intersection of the initial and the final orbit, orbits 
1 and 2. This is feasible if the apogee distance (ra) of one 
orbit is greater than the perigee distance (r,) of the other. 

For the planar problem, y = 0 and Equation [5] becomes 


|AV| = — + — u,*)? (5] 


- The change of velocity |AV | is a function of two variables 
* and which are related by Equation [4]. Therefore, 


AV | is a function of one independent variable, say 3:*, aid 


d(AV)? (AV)? 


In an argument similar to that given in the preceding scc- 


>0 \7 


tion, Loh (7) concluded that the optimum transfer occurs at 


the point of tangency of the two co-planar orbits. The 


argument is not valid in the present case, because when the 
relative orientation changes the location of the point of inter- 
section, changes, and so do the velocities v1, and 


In the following paragraph it will be shown that at the point 
of tangency |AV | is not an extremum. 


At the point of tangency of the two orbits, with J, = 9d, 
and d = + g, the conditions on and are 


+ cos = pill + cos (8 + 


poe: SIN = py €2 Sin + 


| = Que, sin 3; 
dd; Pz 


case of sin 3 = 0, the radial distance of the perigee of one 
orbit is equal to that of the perigee of the other, and 3 = 3; = 


0 or = is the only possible position for the transfer by one 
impulse. 


The value of 3,* is obtained from Equations [7] by Newton’s 
method with # as the starting value. #, does not exist if 
the geometries of the two orbits are such that either they do 
not intersect or they do at all the orientations. For the for- 
mer case transfer by one impulse is impossible, and for the 
latter case the true anomaly of one of the points of intersec- 
tion when 3; = & is employed as the starting value of 3*. 
Numerical data regarding the optimum change of velocity, 
the thrust angle, the location of the orbital transfer and the 
relative orientation of the two orbits are given in the tables 
and the charts of (8). 


Optimum Transfer by Two Impulses 


The two impulses are applied at the point of intersection of 
the transfer trajectory with the initial orbit 1 and the point 
of intersection of the transfer trajectory with the final orbit 
2, respectively. If the specific impulse does not change, t!1e 
requirement of minimum fuel consumption is equivalent to 
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the requirement that the sum of the absolute value of the 
changes of velocity owing to impulse is a minimum, i.e. 


S = |AV,| + |AV2| = min [8] 


where 
AV, = change of velocity from orbit 1 to transfer tra- 
jectory 
AV, = change of velocity from transfer trajectory to orbit 
2 


S is a function of four independent variables, namely, the 
geometry of orbit 3, es and p3and the orientation of orbits 2 and 
3 with respect to orbit 1. The conditions for an extremum 
of % are in general four nonlinear simultaneous algebraic 
equations. In (9), it has been shown that the extremal condi- 
tions are fulfilled by the solutions of the following equations 
(Fig. 2) 


0 


sin 3; = sin 3, = sin 33 = sin 4 [9] 


where 3, and #3 represent the true anomalies of the point of 
intersection of orbits 1 and 3, and 8; and #4 represent the true 
anomalies of the point of intersection of orbits 2 and 3, respec- 
tiv-ly. This statement will be verified in a different manner 
in this paper. It will be shown that the first variation of 
transfer solution subjected to Equation [9] vanishes; the 
second variation is examined to find out whether the extremal 
solition is a local minimum solution. The examination of 
the second variation is simplified when two of the four inde- 
pendent variables therein are uncoupled from the other two. 
This is accomplished by selecting 31, 32, 8; and 3, as the four 
independent variables of S. 

lsefore the study of the extremal solutions, all the possible 
combinations of the geometries of the initial and the final or- 
bit will be classified as follows: 

I] Orbits 1 and 2 never intersect each other 


(ra)2 > (Tp)2 > (%a)i > (Tp)1 
II Orbits 1 and 2 always intersect each other 
(re)i > (Ta)2 > (rp)e > [10b] 


III Orbits 1 and 2 intersect each other at certain orienta- 


[10a] 


tions 
(ra)o > (Tali > > (Tp) [10c¢] 
where 
ra = radial distance of apogee = p/(1 — e) 
ry = radial distance of perigee = p/(1 + e) 


It will be understood that these three combinations also 
include the corresponding ones where the subscripts 1 and 2 
are interchanged. For an optimum transfer trajectory from 
orbit 1 to 2 is also an optimum one for the reverse transfer. 

The combinations of the geometries of the initial and final 
orbits depend on three parameters, p:/ps2, and For a 
given value of é€, the combinations I, II and III lie in the 
corresponding regions in the p;/ps, é: plane as shown in Fig. 3. 

For region I, the condition (rp)2 > (ra), yields the upper 
boundforpr 


2) 


For region II, the lower and upper bounds of p; are 


e P2 + 
It is clear that 


or region III, the lower and upper bounds of p; are” 


l+e pe (1 — e:)/(1 — fore <e 
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REGION I, (fq), > (Tp), >(1g), >(rp), 
REGION IZ, (fq), > (fa), >("p),>('p), 
REGION IZ, (Ta), > (Ta), >("p),>("p ), 


Fig. 3 Classification of combinations of initial and final orbits 


APOGEE-PERIGEE TRANSFER 


— PERIGEE -APOGEE 
TRANSFER 


— APOGEE -APOGEE 
TRANSFER 


7 


6,77 


0 
~~” 
— PERIGEE~PERIGEE TRANSFER 
8,=0, $,=0, ,=0, 


Fig. 4 Combination I—initial and final orbits cannot intersect 
each other 


The solutions of Equation [9] for each of these three com- 
binations are examined separately in the following. 


Combination I—Initial and Final Orbit Never Intersect 
Each Other 


Owing to the restriction of Equations [10a] or [lla], the 
apogee of the transfer orbit 3, corresponding to the solution 
of Equation [9], has to coincide with the apogee or the perigee 
of orbit 2, whereas the perigee of orbit 3 coincides with that 
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of orbit 1. There are four solutions for the transfer from orbit 
1 to orbit 2 (Fig. 4): 

1 Perigee to apogee—orbits 1, 2 and 3 are in same sense. 

2 Apogee to perigee—orbit 3 is in opposite sense to that of 
orbit 1 and 2. 

3 Perigee to perigee—orbit 2 is in opposite sense to that 
of orbit 1 and 3. 

4 Apogee to apogee—orbit 1 is in opposite sense to that 
of orbits 2 and 3. 


For each solution, its first and second variation are examined. 
For the perigee to apogee transfer, the values of es and ps 
are given by the equations 


Ji 1 


The true anomaly of the point of transfer from orbit 1 to 
orbit 3 is zero and that from orbit 3 to orbit 2 is x with g: = 
ge = 0, where ¢, and g, denote the true anomaly of the perigee 
of orbit 3 with respect to orbits 1 and 2, respectively. 

_ The change of velocity at the perigee and the apogee are 


1 1+ 


The small variation of this solution can be represented in 
general as (Fig. 5) a transfer orbit of semilatus rectum p;’ 
and eccentricity es’ with 

[14a] 


ps’ = ps(1 + a) 


intersecting orbit 1 at (7, 6) and orbit 2 at (7,  — 
The true anomalies of the corresponding points with respect 
to the transfer orbit are, respectively, 6; and  — 6. The 
angles 6; — 6, and 6s; — 6. — (6; — 6,) represent the new 
orientations of the transfer orbit and orbit 2 with respect 
to orbit 1. 

With 6;< 1 andj = 1, 2, 3, 4, o and e can be related to 
6; by Equations [1 and 12] 


_1l+e'coss 


e<l 


1 


it 


Fig. 5 Small variation of a solution of Equation [9] 
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1 1+ e:cos (x — , 


1 + cos — 64) 


Ds Ds 2ps Ds 


Therefore 


IIS 


2 


It is evident that ¢ and o are of the order of 6,7. 
At the point of transfer from orbit 1 to the transfer orbit, 
Equations [2] for the radial and circumferential velociti:s 


become 
— - = 
Vu Vn 
(us — (vps! Vp) of, + K, 


Vu ri 
where I, is defined by Equation [13a] and 

As VPs (5) ( Ps pr) [1 
The change of velocity normalized by the factor Vv "i is 

| V [(us — + (vs — 01)?)/u 
+ (LK, + (vs — | 
and the variation in the change of the normalized velocity is 

|| = K,(sg 1) + — [17] 


where (sg I,) means sign of I). 
Similarly, the variation in the change of the normalized 
velocity during the transfer from orbit 3 to orbit 2 is 
I,’ | = IT. | = K, (sg I.) (vs v2)?/( 


vs eyb2/ V ps 


(17b] 


where 


Vs e364 iV ps 


Ke = {(Vp2 — — o(1 — €3)/W/ps} /2 [1s] 


The total variation is ee 
1 f(v2 — v4)? (11 — 
= Kalsg + Kp(se 1 = 
= + 2616163 + + 
21, | 
(254? + + C262”) [19] 
2 |I. | 


In the derivation of Equations [12 to 19], the restriction 
of Equation [10a] has not been imposed; therefore, they are 
valid for any combination of the initial and final orbits. 

The expression on the right-hand side of Equation [19] is a 
linear combination of 6;6;.. It is a second variation since the 
first variation of |I,| + |I,| vanishes. It will be explained 
later how to change the sign of e so that Equations [12 to 19] 
hold for other solutions of Equation [9]. The disappearance 
of terms of first variation verifies that solutions of Equation 
[9] are extremal solutions. 

The uncoupling of 6; and 63; from 64 and 62 obviously simpli- 
fies the study of the second variation. 

For the perigee-apogee transfer between two orbits sul)- 
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p2 > ps > fr, and sg 1, and sg I, are both equal to 1. With 


the definitions of the six coefficients in Equation [19], it 


follows that 
a = (1 é3)(1 + €1) |/ps [20a 


aye, — by? = + e3)(3 — ps V pi)?/(4pips?) 
] 


dy = + e3 + Vil — €3)(1 — é2)]/ps 


ay. — bs? = — €3)(3 + €3:)(W pe — Vps)?/(Apops?) 
[21b] 


It is obvious that a > 0, acy — 6,2 > 0, a2 > O and ac, — 
> for all possible combinations of e:, p: and sub- 
jected to the condition of Equation [lla]. Consequently, 
the second variation of |I, | + I, | given by Equation [19] 
is positive definite, and the perigee-apogee transfer represents 
a local minimum solution. 

ior the apogee-perigee transfer, Equations [12 to 20] hold 
if +; and are replaced respectively by —e, and Con- 
sequently, ajc; — 6:7 is less than 0 and asco — b,? is less than 
0, and the apogee-perigee transfer does not represent a mini- 
mim solution. 

lor the perigee-perigee transfer, all those equations hold 
if ¢ is replaced by —e:. Consequently, aoc. — 0»? is less 
than 0, and the perigee-perigee transfer does not represent a 
minimum solution. 

lor the apogee-apogee transfer, all those equations hold if 
é is replaced by —e,. Consequently, ac: — 6;? is less than 
0, and the apogee-apogee transfer does not represent a 
minimum solution. 

In conclusion, only the perigee to apogee transfer represents 
a local minimum solution for transfers from orbit 1 to 2 which 
do not intersect each other at any orientation. 


Combination Il—Initial and Final Orbit Always Intersect 
Each Other 


Owing to the restriction of Equation [10b], the extreme 
solutions for the transfer from orbit 1 to 2 are perigee to 
apogee, apogee to perigee, perigee to perigee and apogee to 
apogee as shown in Fig. 6. The orientation of the transfer 
orbit for the apogee to perigee transfer and that for the 
apogee to apogee transfer are 180 deg out of phase from the 
corresponding ones in Fig. 4. 

For the perigee to apogee transfer, Equations [12 to 19] 
hold. From Equations [12 and 10b] it is clear that 


[22a] 


After the evaluation of the six coefficients in Equation 
|19], the following is obtained 


a, = 1 + es + V(1 + e5)(1 + [23a] 


aye) — = eyes(1 + €3)(3 — e3)(Wp, — 
[23b ] 


ps < Pr Ds < pe 


sg I, = -1 wi,= 


= V(1 — es)(1 — + (1 [24a] 


ac. — be? = ees(1 — €3)(3 + (Wp, — 
[24b] 


It is clear that a, > 0, az > 0, aye) — b,2 > O and ayes — bo? > 
0, and the perigee-apogee transfer represents a local minimum 
solution. 

For the apogee-perigee transfer, Equations [12 to 19] and 
23 to 24] hold, if negative signs are added to e, e: and e;3. 
Consequently, the apogee-perigee transfer also represents a 
local minimum solution. 

By changing the sign of e in those equations it can be shown 
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jected to the restriction of Equation [10a], it is clear that ; 


—APOGEE- PERIGEE TRANSFER 
4220 
-—-APOGEE-APOGEE TRANSFER 


—PERIGEE-PERIGEE TRANSFER 


6,70 , $,70 
Fig. 6 Combination II—initial and final orbits always intersect 
each other 


(9) that neither the perigee-perigee transfer nor the apogee- 
apogee transfer represents a minimum solution. 

There are two local minimum solutions, the perigee-apogee 
transfer and the apogee-perigee transfer. The total impulse 
or the sum of the absolute values of the changes of velocities 
corresponding to one of those two minimum solutions will be 
compared with that of the other, so as to determine which 
solution of transfer requires less impulse. 

For the perigee to apogee transfer, the total impulse re- 
quired is proportional to I,-. which, according to Equations 
[12, 13 and 22], is 


Vn Vs 
For the apogee to perigee transfer, the geometries of the 
transfer orbit are given by the equations 


Ps P2 ps Pi 


It is evident that py > p; and ps > po, and the total impulse 
required is proportional to I,_, which is 

Pr 
V ps V ps 
The difference of impulse required is proportional to — 

2 2 2 2 a 
Vn Vs Vi 
From Equations [10b, 12 and 26] it can be shown that 


Psps Pip2 


I,-. — = 


PERIGEE -APOGEE a 
TRANSFER 
6,= 7, $220 
7 
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is 
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Therefore, the perigee-apogee transfer is not a minimum 
solution as compared to the apogee-perigee transfer. For 
the latter, the apogee of the transfer orbit coincides with that 
of orbit 1. Since (re); > (ra)2, it is concluded that the trans- 
fer from the apogee of the orbit with longer apogee distance 
to the perigee of the other orbit represents a minimum solu- 
tion. This statement is also applicable for the optimum 
transfer between two orbits of combination I. 


Combination III—Initial and Fina] Orbit Intersect Each 
Other at Certain Orientation 


Owing to the restriction of Equation [llc], p: is less than 
pe. The four extreme solutions are shown in Fig. 7. The 
orientation of the transfer orbit for the apogee-perigee 
transfer is 180 deg out of phase from the corresponding one 
in Fig. 4. 

For the perigee to apogee transfer, Equations [12 to 21] 
hold, since pe > ps3 > pi, and the transfer represents a local 
minimum solution. 

For the apogee to perigee transfer, those equations again 
hold if negative signs are added to &, €2 and €3. This transfer 
also represents a local minimum solution. 

By changing the sign of e in Equations [12 to 19, 23 and 
24], it can be shown that neither the perigee-perigee transfer 
nor the apogee-apogee transfer represents a minimum 
solution. 

Again there are two local minimum solutions, the perigee- 
apogee transfer and the apogee-perigee transfer. The corre- 
sponding transfer orbits are designated as orbits 3 and 4, re- 
spectively. With ps and p, both lying between p, and jn, 
the difference of the impulses required corresponding to these 
two solutions is proportional to 


2e3 
A =1,_. — Il. = — + = = 
A is a function of @, 2, p: and pe. Proving that it is positive 
or negative in a straightforward procedure similar to Equa- 
tions [27 and 28] has not been. successful. Instead, it is 
shown in (9) that 0A/dp: > 0. Along the upper bound of 
71, according to Equation [llc], either the perigee or the 
apogee of orbit 1 coincides with that of orbit 2 and A vanishes. 
Therefore in region III of the p;/pe, é2 plane (Fig. 3) A < 0; 
i.e., the apogee-perigee transfer is not a minimum solution as 
compared to the perigee-apogee transfer. For the latter, 
the apogee of the transfer orbit coincides with that of orbit 2. 
Since (ra)2 > (re): according to Equation [10c], the general 
statement regarding the optimum transfer between two orbits 
of combination II applies also for the present combination. 
Concluding Remarks 


It has been shown that the minimum impulse transfer from 
one given orbit to another of given geometry can be accom- 
plished only when the two orbits and the transfer trajectory 
are co-planar. 

If the two orbits can intersect each other at some relative 
orientation, the orbital transfer can be executed by a single 
impulse. The procedure of finding the optimum orientation 
has been presented. 

The transfer from a given elliptical orbit to another of given 
geometry by two impulses has been studied. It is shown 
that the total impulse will be a minimum, provided that the 
two orbits and the transfer trajectory are co-planar, with 
the same major axis and oriented in the same sense, and 
provided that the apogee of the transfer trajectory coincides 
with that of the orbit, with longer apogee distance, and the 
perigee of the transfer trajectory coincides with that of the 
other orbit. This optimum transfer trajectory represents a 
generalization of Hohmann’s transfer between two concentric 
circular orbits. 

It should be pointed out that the solutions of Equation 
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Fig. 7 Combination [[]—initial and final orbits can intersect 
each other at some orientations 


[9] are extremal solutions, but they may not necessarily repre- 
sent all the extremal solutions. Investigation in this respect 
will be carried out. 

If the two orbits can intersect each other, orbital transfer 
can be accomplished by either one or two impulses. By 
comparison of numerical results, the sum of the two impulses 
inthe optimum transfer is always less than the corresponding 
optimum single impulse. An analytic proof of this statement 
is desirable. Although the transfer from one orbit to another 
can always be accomplished by two impulses, Shternfeld (10) 
and Hoelker and Silber (11) have shown that the transfer 
between two circular orbits can be accomplished more 
economically by three impulses if the ratio of the larger 
radius to the smaller one is greater than 11.9. As a generaliz- 
ation of their investigation, the optimum transfer between 
two elliptical orbits by three impulses or more has been 
studied (12). 
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numbers above approximately 


N COMPARING the performance of conventional ram- 
jets and ramjets having supersonic combustion, the Mach 
number of combustion and the total pressure losses in the 
inlet diffuser for the supersonic combustion ramjet should be 
identical to the Mach number upstream of the inlet normal 
shock and the total pressure losses owing to diffusion to this 
Mach number, respectively, for the subsonic combustion 
ramjet. Although several comparisons of supersonic combus- 
tion and subsonic combustion ramjets were made in (1),? 
none were made for the case where each ramjet type employed 
inlets having the same minimum supersonic Mach number 
and the same total pressure losses in the supersonic diffusion 
process. 

The purpose of the study reported here was to determine 
the relative performance of the two ramjet types for the two 
cases where each system employs an isentropic spike external 
compression inlet and where each system employs an internal 
compression inlet. The study is based on a simple analytical 
model rather than on complete cycle studies such as those 
used in (1). 


Discussion 


Typical Mach numbers which might exist in a supersonic 
combustion hypersonic ramjet are shown in the upper part of 
Fig. 1. The Mach number at the inlet throat (14) was 
chosen so that the addition of fuel at an equivalence ratio 
¢ of 0.5 would result in sonic flow at the burner exit. 

Also shown in Fig. 1 is a subsonic combustion hypersonic 
ramjet operating with the same inlet throat Mach number, 
the same free stream Mach number (Mo) and the same equiv- 

alence ratio as the supersonic combustion ramjet shown at 
the top of the figure. The geometry of the subsonic combus- 
tion ramjet is unconventional in that the pressure rise through 
the normal shock and the combustion of the fuel take place in 
the same constant-area passage. If wall friction in the 
constant-area section is neglected, then the same input 
information (pressure, temperature and _ velocity) and 
equations (continuity, momentum and energy) are employed 
to calculate the change of properties across the combustor in 
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Relative Thermodynamic Efficiency of 
Supersonic Combustion and Subsonic 
Combustion Hypersonic Ramjets 


An analytical model is set up to permit comparison of the relative thermodynamic efficiencies 
(fuel consumption or specific impulse) of supersonic combustion and subsonic combustion hyper- 


sonic ramjets for arbitrary inlet performance. 
supersonic combustion ramjet will be superior to that of a subsonic combustion ramjet at Mach 
7 for configurations having isentropic spike external compression 


inlets and above 10 for configurations having internal compression inlets. .- 5 cee 


and subsonic combustion hypersonic ramjets. 
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The results indicate that the specific impulse of a 


4 


the supersonic combustion ramjet as across the normal shock 
and combustor of the subsonic combustion ramjet. Usually 
there are two solutions to these equations, one indicating 
supersonic flow at Station 3 and the other indicating sub- 
sonic flow at Station 3. However, since the equivalence ratio 
was chosen to give a Mach number of unity at Station 3, 
the two solutions are identical, the conditions at Station 3 
(pressure, temperature and velocity) are the same for both 
ramjets, and both ramjets would produce the same thrust with 
the same specific impulse or fuel consumption. 

The effect of variations of equivalence ratio from the value 
employed in Fig. 1 on the ratio of burner exit total pressure to 
free stream total pressure (which determines specific im- 
pulse) for a fixed value of M, is shown in the upper part of 
Fig. 2. A decrease in equivalence ratio results in a relatively 
small reduction in the total pressure losses in the subsonic 
combustion ramjet, since the majority of the total pressure 
losses occur across the normal shock wave for this configura- 
tion. (A reduction in equivalence ratio requires a reduction 
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in area at the burner exit to maintain proper matching of 


: components in the subsonic combustion ramjet.) However, 


_ for the supersonic combustion ramjet, a reduction in the 
heat addition in the supersonic stream results in a large 
decrease in total pressure losses across the supersonic 
combustor and produces a supersonic rather than a sonic 
velocity at the burner exit. Conversely, the increase in 


My = 8 


SUPERSONIC-COMBUSTION RAMJET 


SUBSONIC- COMBUSTION RAMJET 


EQUIVALENCE RATIO, 


TOTAL PRESSURE AT BURNER EXIT 


INLET THROAT OR BURNER ENTRANCE MACH NUMBER,M, 


Fig. 2 Effect of equivalence ratio and inlet throat Mach number 
on total ames losses in supersonic combustion and subsonic 


combustion hypersonic ramjets 
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total pressure losses resulting from an increase in equivalence 
ratio is much less for the subsonic combustion ramjet than 
for the supersonic combustion ramjet. If the equivalence 
ratio is increased above 0.5. it is necessary as a result of 
continuity considerations to have the combustion from the 
added fuel take place in an increased cross-sectional arca. 
In the supersonic combustion ramjet, this additional coin- 
bustion would take place in the divergent portion of thie 
nozzle, thereby creating additional losses in total pressure. 
For the subsonic combustion case, however, it would be more 
advantageous to increase the area of the combustion ¢ _ er 
between the inlet throat and the exhaust nozzle throat } 

order to allow combustion of the increased amount of fuel. 

Although the reasoning presented in the preceding paia- 
graph was given in terms of a fixed Mach number at the in ct 
throat and a variable equivalence ratio, the same kind of 
reasoning can be used for a fixed equivalence ratio and a 

variable inlet throat Mach number (see lower half of Fig. -). 
If the Mach number (M;) provided by the inlet is greater 
than that which would result in M; = 1, then the supersonic 
combustion configuration is more efficient. Conversely, 
the Mach number provided by the inlet is less than that which 
would result in M; = 1, then the subsonic combustion con- 
figuration is more efficient. 

Since the breakeven free stream Mach number for super 
sonie combustion and subsonic combustion ramjets is a 
function of the Mach number delivered by the inlet to tiie 
throat of the configuration, it is of interest to examine tlie 
throat Mach numbers which would exist for various types of 
supersonic inlets. It can be shown that it is always desirable 
for both subsonic combustion and supersonic combustion 
Ffamjets to decrease the inlet throat Mach number as much «as 
possible. (A decrease in M, in a subsonic combustion ram- 
jet decreases the total pressure losses across the normal 
shock. A decrease in M, in a supersonic combustion ramjct 
decreases the total pressure losses in the combustion process, 
even though a decrease in M, in some cases might require that 
some of the combustion take place downstream of Station 3.) 
However, the minimum attainable value of M, for supersonic 
inlets is limited by various aerodynamic considerations. ‘The 
variation of minimum attainable inlet throat Mach number 
with free stream Mach number is shown in Fig. 3 for three 
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Fig. 3 Comparison of minimum attainable inlet throat Mach numbers with critical burner entrance Mach numbers 
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types of supersonic inlets: External compression inlets with 

no internal contraction, external compression inlets with 
internal contraction, and internal compression inlets. The 
curve shown for the external compression inlets with no 
internal contraction was obtained from a procedure which is 
similar to that outlined in (2), but which accounts for real 
gas effects and boundary layer growth. The minimum 
attainable Mach numbers for the external compression inlet 
with internal contraction shown in Fig. 3 were obtained by 
reducing those shown for the inlet with no internal contrac- 
tion by the maximum allowable amount for fixed-geometry 
inlets. The range of minimum attainable throat Mach 
numbers shown for internal compression inlets was deter- 
mined on the basis of experience obtained with these inlet 
types at the UAC Research Laboratories. 

Also shown in Fig. 3 are critical values of burner entrance 
Mach number (i.e., values of burner entrance Mach number 
M, for which the addition of heat results in a Mach number of 
unity at the downstream end of the burner). The critical 
values of burner entrance Mach number shown in Fig. 3 
wer’ determined by a procedure which accounts for all real 
gas effects. As noted in a preceding paragraph, the Mach 
number at which supersonic combustion and _ subsonic 
combustion ramjets produce the same thermodynamic effi- 
cienvies is determined by the free stream Mach number at 
which the critical value of burner entrance Mach number 
is equal to the throat Mach number provided by the inlet. 
It is shown in Fig. 3 that the breakeven free stream Mach 
nuniber varies from approximately 7 for external compression 
inlets to approximately 10 for internal compression inlets. 

The critical value of burner entrance Mach number is a 
function of the fuel-air ratio and method of fuel injection in 
the burner. If the equivalence ratio is decreased below 1.0, 
the critical inlet throat Mach number will decrease slightly 
(this decrease is very small at high Mach numbers because the 
fuel added at high equivalence ratios results in an increase in 
dissociation rather than an increase in temperature). The 
critical burner entrance Mach numbers for equivalence ratios 
greater than 1.0 are also slightly less than those shown on 
Fig. 3 for an equivalence ratio equal to 1.0. Upstream 
injection of the fuel results in an increase in critical burner 
entrance Mach number from the values shown in Fig. 3. 

The efficiency of a subsonic combustion configuration could 
be improved over that obtainable with the subsonic combus- 
tion configuration shown in Fig. 1 by increasing the area of 
the passage downstream of the inlet throat and by injecting 
and burning the fuel at the resulting reduced Mach number. 
Therefore, the breakeven free stream Mach numbers for this 
alternate, more conventional, subsonic combustion configura- 
tion would be higher than those given in Fig. 3. 

It has been suggested that the weight of a supersonic 
combustion configuration would be less than that of a sub- 
sonic combustion configuration. Weight is determined by 
three factors: Pressure loads, length and cooling require- 
ments. Although the pressures within a supersonic combus- 
tion configuration are less than those in a subsonic combustion 


configuration, it might be necessary to design the supersonic 
combustion configuration to provide for the possibility that 
a normal shock would be created inadvertently in the inlet 
throat. If this condition were considered, there would be no 
difference in weight owing to differences in pressure loads 
within the two configurations. 

The overall length of a burner at high free stream Mach 
numbers is determined primarily by the length of the fuel 
injection portion of the burner. It seems doubtful that the 
fuel can be injected in a supersonic stream in a shorter distance 
than in a subsonic stream. 

The rate of heat transfer to the ramjet duct wall per unit 
area at a given free stream Mach number is a function of the 
mass flow per unit duct area and the wall friction coefficient. 
Although there is no difference in the mass flow per unit area 
for the two configurations shown in Fig. 1, the friction co- 
efficients in supersonic flow are slightly less than those in 
subsonic flow. Therefore, the heating load in*the super- 
sonic combustion configuration would be somewhat less 
than that in the subsonic combustion configuration if the 
supersonic burner could be made as short as the subsonic 
burner. 

Dissociation of the products of combustion in high Mach 
number ramjets results in a loss in thrust if the reaction rates 
in the exhaust nozzle are not fast enough to maintain chemical 
equilibrium. For the equal-thrust configurations shown in 
Fig. 1, the static temperatures at Station 3 for chemical 
equilibrium at Station 3 are equal for the supersonic combus- 
tion ramjet and the subsonic combustion ramjet. For this 
case the chemical reaction rate in the exhaust nozzle has no 
influence on the breakeven free stream Mach number between 
supersonic combustion and subsonic combustion ramjets. 
If the flow at Station 3 is not in chemical equilibrium, the 
dissociation will be greater for the subsonic combustion ram- 
jet than for the supersonic combustion ramjet, as a result of 
the higher static temperatures upstream of Station 3 for the 
subsonic combustion ramjet. This increased dissociation 
will cause a loss in thrust and serve to decrease the break- 
even free stream Mach number from that determined by the 


construction given in Fig. 3. a 

Nomenclature 

My, = free stream Mach number 

M, = inlet throat or burner entrance Mach number 

Mz, = Mach number behind normal shock in subsonic combus-_ 
tion configuration 

M; = burner exit Mach number 

@ = equivalence ratio, ratio of actual fuel flow to fuel flow for 


stoichiometric mixture 
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fect if Losses on the 
Desi of Absorhers for Combustion 


hy one 5 ee, She Investigations were conducted to provide a basis for the design of multiple-resonator absorber- to 
4 aaee eats pls: control high amplitude combustion instabilities (screech) in high output combustors. An imped- 
fom cane = ane ance tube technique was employed to determine the effect of high amplitude sound pressure waves ) 
oh oi the absorption characteristics of absorber configurations having open areas varying between ‘).5 | 
+ 4 rest i ne and 10 percent. Frequencies were varied between 200 and 1200 cps. It was found that the acousti- | 
aa ah cal impedance of the absorbers contains a nonlinear facing resistance which depends on the oscillat- | 
ee nate ing velocity and particle displacement in the absorber neck, and that the nonlinear resistance can | 


generally be correlated with the amplitude of the incident pressure wave plotted on a decibel scile. 
oe * pe bi nas A design procedure is presented in which the absorption characteristics of perforated liners for a 
hig tess ae particular sound wave amplitude can be calculated from known physical dimensions. 


CREECH is a type of combustion instability encountered Theoretical Considerations 7 
in high output combustors, such as those used in rockets, 
ramjets and turbojet afterburners. It is characterized by Resonator Array 


high amplitude pressure and flame front oscillations which 
occur at one of the resonant frequencies of the combustion 
chamber. These oscillations increase heat transfer and 
stress repetition rates, and usually cause failure of structural 
components of the combustion chamber in a matter of sec- ( 
onds. It is obvious that combustion instabilities of this ( 
type must be avoided. § 
Previous investigations have shown that screech pressure a = 46/[(6 + 1)? + xX] (1] 


Fig. 1 illustrates an array of closely spaced resonators 
typical of configurations used in screech absorption liners. 
If the perforated area is assumed to have a uniform imped- 
ance, the absorption coefficient (defined as the fractional part 
of the energy of an incident sound wave absorbed by the 
liner) can be written for normal incident waves (2) as 


oscillations can be suppressed by lining the combustion : 
chamber with a perforated metal sheet to absorb the kinetic where @ is the specific acoustic resistance ratio and x is the 
energy of the oscillations in much the same manner that a specific acoustic reactance ratio, which can be expressed as 
Helmholtz resonator absorbs sound. The problem of screech Qafol :-% \s i 
liner design is one of adapting the simple Helmholtz resonator ( ( [2] 
theory to the conditions of high pressure amplitude that an fo ff 

exist when screech occurs. Previous investigations of single Hence en 7 

resonator configurations (1)? have shown that absorption is a L 

function of pressure wave amplitude when the amplitude = 46 
level is greater than about 100 db. This relationship is a nas Zrfolere f fo\? [3] f 
obtained at the higher amplitudes because a resonator con- (1 + 8) in y ie -*) | 
tributes a nonlinear resistance to the sound wave in addition 

to the linear resistance. This increase in resistance results which for 6 = land f = fogives a = 1. 
from changes in the circulation and turbulence of the gas 

particles at the neck of the resonator as the pressure wave 7 

amplitude or particle velocity increases. At levels greater 

than about 100 db, the nonlinear resistance becomes a con- a PERFORATED I 
trolling factor in determining the absorption characteristics pats ak 
of aresonator. This resistance has been measured for certain T T 2,1 
limited conditions, but has not been quantitatively expressed T 


co 


for the perforated plate, i.e., the resonator array. a 
The purpose of the investigation described in this paper 7 
was to determine the effects of high amplitude pressure waves a re 
on the absorption characteristics of multiple resonator 
J 


absorbers and to incorporate these results in the development 
of a design procedure for screech absorption liners. 
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aA, Fig. 1 Resonator array 
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Fig. 2a Aperture phase dia- 

gram showing regions of 

acoustic streaming and 


turbulence particle velocity 
vs. plate thickness (1) 
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The assumption of uniform impedance applies only if the 
distance a between the centers of the resonator necks spaced 
in a square array is less than \/(27)!/* (approximately one 
half wave length) and greater than the diameter d of the 
aperture (1,3). For values of a > d/(2m)!/2, the individual 
resonators should act independently, and the total absorption 
should be determined as the sum of the individual absorption 
cross sections of the resonators. For values of a approaching 
d, the resonators should probably be treated as a pure re- 
sistive load with no reactance. The dimensions of all prac- 
tical liners for screech absorption are expected to fall within 
the range where the assumption of uniform impedance 
applies. 

If the pressure waves striking the absorber are not of normal 
incidence (4) 


46 cos 
(6 cos + 1)? + (2s - cos *¢ 


Partitions located in the resonator volume behind the per- 
forated plates are generally necessary if the pressure waves 
are not of normal incidence; however, for normal incident 
waves, partitions have little or no effect on the absorption 
characteristics. 


Effect of High Intensities on Absorption Characteristics of 
Resonators 


Previous investigators have shown that the absorbing 
characteristics of resonators are materially different when 
exposed to high amplitude sound waves rather than low 
amplitude sound waves. The most thorough investigations 
of this phenomenon, as reported in (1 and 5), have shown that 
the flow of air particles in and out of the resonator neck can be 
classified in one of four regions. Phase diagrams of these 
regions are presented in Figs. 2a and 2b. Region 1 is a low 
intensity sound region with stationary circulation; the flow 
is directed out of the orifice along the orifice axis. Region 2 
is a region of stationary circulation in which the direction of 
flow along the axis is toward the orifice. Region 3 is charac- 
terized by pulsing effects superposed on circulation of the 
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kind in Region 2. Region 4 is a high intensity sound region 
in which pulsatory effects predominate and result in the 
formation of jets and vortex rings. 

The losses contributed by the turbulence and circulation 
associated with these flow regions are nonlinear and similar 
to those caused by the facing and aperture resistance; they 
may be incorporated in the equation for resonator resistance 
by a term A,, to give 


( t 5 


where the terms 1, t/d and A,./d represent the facing, neck 
and circulation losses respectively. 

At incident wave pressure levels of the order of 100 db, the 
nonlinear losses represent a significant part of the total 
resistance. When the levels are as large as 130 db, the total 
resistance is almost entirely a result of this term. Since 
screech pressure levels are of such order of magnitude, it is 
essential that the A,: term be accurately known for the design 
of screech liners. Some evidence also indicates that addi- 
tional nonlinear terms should be included if there is flow past 
or through the resonator face (6). In the work described 
here, these losses have been assumed to be small relative to 
the intensity dependent loss and have not been included. 

Some work has been done to determine the nonlinear losses 
as a function of operating variables (1), and the results have 
been expressed in accordance with the phase diagrams previ- 
ously mentioned; i.e., the nonlinear resistance is correlated 
with a variable whose relation to the region boundaries is 
known. It was found that the onset of the turbulent regions 
in the phase diagrams (Figs. 2a and 2b) was a function of fre- 
quency f, particle velocity u and resonator aperture thickness 
t. A critical thickness was defined as 


= u,/2 

te = [6] 
where u- is the velocity corresponding to the lower boundary 
of Region 3. The dashed lines in Figs. 2a and 2b correspond 
to u = 2x ft in which the particle displacement in the aperture 
equals the aperture thickness ¢t. For sufficiently large values 
of t (t > t.), the turbulence onset occurs approximately at a 
constant value of u and is independent of the aperture thick- 
ness t. Ift < t., the criterion for the onset of turbulence is s/t 
= 1, where s is the particle displacement amplitude. The 
critical thickness t. was found to be of the order of 0.0236 in. 
(1). 

The nonlinear resistance was found to be consistent with 
the observations on turbulence. For t < t., the nonlinear 
resistance depends essentially on s/t, but for t > t. the particle 
velocity amplitude is the governing variable. Ingard’s 
results (1) were of limited value when applied to screech liner 
design because multiple resonator effects, frequency variation 
and reactive loading were not considered. It was a purpose 
of the work described here to consider these effects. 


Apparatus and Procedure for Experimental 
Investigation 
Method of Measuring Absorption Coefficients 


The impedance tube technique was utilized to measure the 
absorption coefficients of the various test configurations. It 
can be shown (7) that the absorption coefficient can be deter- 
mined by placing the test configuration at the end of a tube 
which is excited at a resonant frequency by some driving 
source and by measuring the magnitude of the pressure 
maxima and minima of the standing wave pattern in the tube 
and the distance between the pressure maxima and minima. 


Determination of Natural Frequency 


The natural frequency of a resonator configuration is the 
frequency at which the resonator offers no specific reactance 
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Fig. 3 Instrumentation diagram for impedance tube method of 
sound absorption measurement 


xX to the sound wave. This frequency, called the resonant 
frequency, can be expressed as follows for the resonator 
array 

c 


fo = 2r La 


Because the ratio of the minimum pressure amplitude to 
the maximum pressure amplitude of the standing wave in the 
impedance tube is a maximum at resonance, it was deter- 
mined that a plot of pmin/Pmax vs. frequency was a suitable 
method for obtaining a value of fy from experimental data. 
A plot of the minimum pressure position vs. frequency 
was used as a check. 


All configurations were designed so both the backing dis- 
tance and the position of the test configuration in the imped- 
ance tube could be varied. Aperture plates were cut from 
z¢-in. steel sheets. The per cent open area of the test con- 
figurations varied between 0.5 and 10 per cent. Apertures 
were spaced in a symmetrical square array so the openings 
had an equivalent absorbing or working area. Apertures 
could not be spaced around the outer portion of the plate 
because of the lack of symmetry which this would have 
produced. However, additional holes were added around 
the edges of one configuration to determine if these holes 
changed the absorption characteristics appreciably. There 
did not appear to be a significant difference between the two 
cases. 


Impedance Tube 


The impedance tube was a 4-in. diameter brass tube about 
4 ft long. One end of the tube was closed with the perforated 
plate absorption configuration. At the other end, a loud- 
speaker was placed which propagated sound waves of con- 
trolled frequency and intensity. The incident wave added 
to the sound wave reflected from the surface of the absorber 
to result in a standing wave along the tube. A traveling 
microphone inserted in a slit cut along the upper length of the 
impedance tube indicated the position and amplitude of the 
maximum and minimum pressure points from which the 
absorption coefficient was determined. 
The resistance of the tube was measured and w as s found 
to cause negligible attenuation. 


The test equipment used for this investigation is illustrated 
in Fig. 3. The frequency was determined by means of a 
Berkeley EPUT meter and a Hewlett-Packard oscillator. 
For a particular sound pressure level intensity the frequency 
was set so the pressure amplitude at the plate would be 
maximized. The maximum pressure amplitude of the 


standing wave was adjusted by means of a Fisher amplifier 
and measured with a General Radio condenser microphone 
pickup system and a Hewlett-Packard harmonic wive 
analyzer. A Ballantine voltmeter was used to indicate the 
oecurrence of a maximum pressure amplitude at the plate. 
The voltmeter was not used for recording actual pressures 
because it lacked the filter medium of the analyzer. Micro- 
phone pickup tubes of different diameters were used in or: er 
to measure pressures accurately over a large range of :n- 
tensities. 


Discussion of Experimental Results 


The absorption characteristics of test configurations having 
open areas varying between 0.5 and 10 per cent were det:r- 
mined for various frequencies at pressure levels between | 00 
and 160 db. The impedance tube technique was utiliz d; 
measurements allowed computation of the real and imagin:.ry 
parts of the impedance for the various test conditions. 


Effect of Intensity Level on Absorption at Resonance 


From experimentally determined values of the reson:nt 
frequency and facing resistance, it was possible to compiite 
the nonlinear resistance factor A,./d at the various pressiire 
levels from Equation [5]. Following Ingard’s procedure ||), 
attempts were made to correlate the nonlinear resistance 
factor with both the particle velocity and particle displace- 
ment in the neck of the resonator. The particle velocities 
and displaceme nts were calculated from the following equa- 
tion (2) 


2P; 
u = 2xfs = |- 

c}(@+1 

Fig. 4a shows that the nonlinear resistance can be correlated 
with the particle velocity of the liner configuration. ‘lhe 


method of least squares was used to determine the location 
of the solid line through the experimental data; the dashed 
line corresponds to the data obtained for a single resonator 
configuration (1). The solid line can be expressed empirically 
as 


logio (An/d) = —0.009 + 0.400 logio (u) [9] 


It can be seen that some divergence exists between the slope 
of the two lines. Presumably this results from the coupling 
of the multiple resonators. The particle velocities were 
usually relatively high, and the acoustic jet (Region 4 of 
Figs. 2a and 2b) should have been formed under most condi- 
tions. Ingard (1) indicates that the nonlinear resistance 
should correlate with particle velocity if the aperture thick- 
ness is above a critical value and with the ratio of the particle 
displacement to the aperture thickness if the aperture thick- 
ness is below the critical value. Although the aperture thick- 
ness used in this investigation was greater than the critical 
thickness given by Ingard (1), Fig. 4b shows that some degree 
of correlation does exist between the nonlinear resistance and 
the particle displacement. Again, the method of least 
squares was used to determine the location of the solid line 
through the experimental data points; the dashed line 
corresponds to the data given by Ingard (1) for a single 
resonator configuration. The solid line is given by the 
empirical equation 


logio )- 0.516 + 0.482 logis (10] 

The scattering of points indicates that the nonlinear re- 
sistance is really a complex function of many variables, but 
it will be shown later that the degree of correlation which 
exists is adequate for design purposes. 
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By plotting the sound pressure level vs. the nonlinear re- 
sistance, a fair degree of correlation was obtained as shown 
in Fig. 5. The solid line through the data points was deter- 
mined by the method of least squares and can be expressed as 


logio (An/d) = —1.685 + 0.0185 (db) [11] 


The degree of fit is somewhat misleading, since the pressure 
amplitude doubles with a mere 6-db increase. However, 
this presentation is a satisfactory and convenient correlation 
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RATIO OF PARTICLE DISPLACEMENT AMPLITUDE IN THE APERTURE TO APERTURE THICKNESS , s/t 


Fig. 4b Nonlinear resistance vs. particle displacement am- 
plitude in the aperture 


for liner design purposes because the absorption coefficient is 
rather insensitive to variations in the nonlinear resistance 
term. For the accuracies presently required for screech 
absorber design, it is suggested that the correlation of Fig. 5 
be used for determining the nonlinear resistances for per cent 
open areas between approximately 0.5 and 10 per cent. The 
use of Fig. 5 for design purposes is discussed in more detail 
in the subsection titled “Improvements in Accuracy of 
Design System.” 


Effect of Intensity Level on Absorption at Nonresonant 
Conditions 


The correlations presented previously apply rigorously only 
at resonant conditions where the reactive part of the imped- 
ance is zero and the mass in the aperture neck vibrates in 
phase with the incident wave. It is of some practical interest 
to be able to determine the absorption characteristics of 
screech liners at off-resonance conditions, for the liner natural 
frequency does not necessarily correspond to the screech 
frequency for an optimum liner design. The reactance of a 
multiresonator absorber can be computed for various fre- 
quencies from Equation [2], and the absorber facing resistance 
can be computed from Equation [5]. The absorption 
coefficient can be expressed in terms of the reactance and 
resistance using Equation [1]. Figs. 6a and 6b_ present 
experimental data on specific acoustical reactance and 
resistance typical of data obtained for the various liner 
configurations. Fig. 6a shows that the reactance increases 
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with frequency and decreases slightly with increasing sound 
pressure level. This result is in agreement with Ingard’s 
findings (1). 

Fig. 6b shows the variation of the facing resistance 6 with 
frequency. From Equation [5] it can be seen that if A,./d 
were independent of frequency, @ would vary with f!/2._ This 
is evidently not the case. No consistent relationship could 
be found between the facing resistance 6 and frequency for 
the various test configurations at the high sound pressure 
levels. In some instances 6 increased with frequency and in 
others it decreased. At frequencies far from resonance, the 
resistance became small. 


Effective Length of Vibrating Air Mass in Resonator Neck 


In the computation of the natural frequency and reactance 
of resonator configurations, it is necessary to know the effec- 
tive length of the vibrating mass. An end correction applied 
to the aperture neck makes the length of the vibrating air 
mass in the neck greater than the plate thickness ¢ (sec 
Fig. 1). As discussed by Ingard (1), the effective length can 
be written as 


lete = t + 0.85d(1 — 0.7 Vc) — 6 [12] 


where 6 is an experimentally determined correction factor. 
Ingard (1) found that 6 could be correlated with the sound 
pressure level for single resonators. 

Values of the correction factor were determined by means 
of Equations [7 and 12] from measurements of the natural 
frequency of the resonator configuration for sound pressure 
levels between 100 and 160 db. It was found that 6 could 
not be correlated directly with the sound pressure level, but 
the per cent open area of the configuration had a predominant 
effect on 6 as shown in Fig. 7. The value of 6 appears to 
decrease until open areas of around 6 per cent are reached 
and then remains essentially constant at higher per cent open 
areas. 


Design System for Screech Absorption Liners 


In presenting a design system for screech liners, some con- 
sideration must be given to the required accuracy and utility 
of the system. With this in mind, a basic system is presented 
in which accuracy is compromised slightly for the sake of 
utility, and methods for improving the accuracy of the system 
at the cost of increasing the tediousness of the problem are 
indicated. 


Basic Design System 


An examination of the experimentally determined curves of 
acoustical reactance and resistance vs. frequency (see Figs. 6a 
and 6b) shows that a reasonable approximation to the curves 
can be made by assuming that the reactance is independent 
of the sound pressure level, and assuming that the resistance 
varies with sound pressure level but is invariant with fre- 
quency. On the basis of these assumptions, the equations 
describing the operation of the absorber can be written as 


46, 


+ 1)? + (F/fo — fo/f)? 
where 
Q = 2nfolest/cob [14] 
4(upfo/g) /2(] + t/d + 
[15] 
apc 


0 

less = t + 0.85d(1 — 0.7 Vo) — 6 
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In most practical liner designs the backing distance L and the 
liner thickness ¢ will generally be fixed because of consider- 
ations of strength and cooling requirements. The diameter 
of an aperture that can be conveniently punched in a plate of 
thickness ¢t, as well as cooling flow requirements, put limita- 
tions on the values of d which can be tolerated. With these 
restrictions, the per cent open area o appears to be the most 
convenient variable to control in designing a liner. 

Equations [13 through 17] can be solved for a selected 
value of by using Figs. 5 and 7. An optimum can then be 
determined by plotting the absorption coefficient vs. selected 
values of frequency f with o as a parameter. The optimum ¢ 
is that value which gives the maximum value of the absorp- 
tion coefficient over the widest frequency bandwidth. The 
same procedure can be followed to optimize the aperture size d 
if there is sufficient variation allowable in d to warrant the 
added difficulty. 

In cases where the absorption bandwidth is too narrow, it 
will be necessary to change the value of Q (see Eq. [14]) for 
the system. The Q value is proportional to the inverse of 
the damping coefficient and is a measure of the absorption 
bandwidth. A system having a high Q value absorbs very 
selectively and has a narrow bandwidth. Systems with a low 
Q value usually have less peak absorption but absorb over a 
wider range of frequencies. Therefore, judgment must |e 
used to achieve a required absorption over the required frv- 
quency range. Variations in Q can probably best be achieved 
by varying the aperture diameter d, and both o and d should 
be optimized if the required bandwidth cannot be achieved by 
optimizing alone. 


ARS JouRNAL 


~ 
~ 
baie. 
> 
] 
Bs 
- 
-~ | r 
a 
~ ‘ 
ti 
{ ul 
at 
In 
al 
lin 
. 
lit 
In 


= 

d= 0.0768 in. 

\ 

=1.45% 

08 T=1.45% 
a 
u 
z 
04 
a 
a 


/ 
EXPERIMENTAL = 


i 
a2 DETERMINATION 
128 db 
/ 
‘4 


CONSTANT 
VARIABLE 


(0) 200 400 600 800 1000 1200 
FREQUENCY, f - cps 


Fig. 9a Comparison of measured and calculated absorption 
coefficients—128 db 


10 T 
t= 0.0625 in. 
d=0.0768 in. 
L=1.9 in 
as T=1.45% / 
‘a 
u Q6 
/ 
/ 
2 oa 
fe) / 
/ 
a © EXPERIMENTAL 
0.2 L 4 DETERMINATION 
\/o 148 db 
@ CONSTANT 
VARIABLE 
fe) 
Ce) 200 400 600 800 1000 1200 


FREQUENCY, f - cps 


Fig. 9b Comparison of measured and calculated absorption 
coefficients—148 db 


If variations in d are impractical, the effective Q of the 
system can be controlled by providing multiple absorption- 
liner sections having different peak frequencies. 

The simplifying assumptions introduced in the design sys- 
tem appear to be justified because large errors in the value of 6 
result in small errors in the absorption coefficient. In Fig. 8 
a comparison is presented of the calculated and experimental 
values of the resonant absorption coefficient as a function of 
sound pressure level obtained from a typical test configura- 
tion. Good agreement was obtained, and the effect of the 
nonlinear losses is quite apparent. The nonlinearities are 
unimportant below sound pressure levels of about 80 db, but 
at higher levels they change the absorption coefficient mark- 
edly. At approximately 140 db the nonlinear corrections 
increase the value of 6 to 1, and a theoretical 100 per cent 
absorption is obtained. Further increase in 6 by the non- 
linear corrections causes the absorption coefficient to de- 


crease. The importance of the nonlinear considerations in 
liner design is obvious. altel 
Improvements in Accuracy of Design System 


Sources of error in the basic design system are the assump- 
tion that the facing resistance is independent of frequency 
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and the use of the correlation presented in Fig. 5 to determine 
the nonlinear resistance term A,,/d. Without these simplify- 
ing assumptions, Equation [15] in the equations describing 
the operating characteristics of the absorber is replaced by 
Equation [5]. (It should be noted that f and not fy is used in 
Eq. [5] for computing 6.) Equation [17] can be solved for 
lere as before by using Fig. 7, and Equation [16] can be 
solved for fo. If the correlations of Figs. 4a and 4b are used, 
it is necessary to calculate the particle velocity and/or the 
displacement in the resonator neck. The particle velocity 
can be computed from Equation [8], where P; is the pressure 
amplitude of the incident wave. The particle displacement 
is given by ~ 


s = u/2nf [18] 


The correlation of Fig. 4a should be used to determine 
A,./d except for very thin plates where the plate thickness is 
less than the critical value of approximately 0.023 in. (1). 
For these thin plates, the correlations of Fig. 4b should be 
used. For plate thicknesses greater than critical, a given a, 
sound pressure level, aperture diameter and frequency, Equa- 
tions [5 and 8] can be solved by trial and error by using Fig. 
4a. The absorption coefficient can then be determined from 
Equation [3]. The process can then be repeated for various 
frequencies to determine the absorption coefficient as a 
function of frequency. If the plate thickness is less than 
critical, the correlation of Fig. 4b can be used in a similar 
manner in conjunction with Equations [5 and 18]. 

Figs. 9a and 9b present a comparison of experimentally 
measured absorption coefficients obtained for a typical con- 
figuration and absorption coefficients calculated using both 
the simplified (6 constant) and rigorous (6 variable) pro- 
cedures. Good agreement was obtained between the meas- 
ured and calculated values. The more complicated calcula- 
tion procedure described here appears to give better agree- 
ment with the experimental data at the higher sound pressure 
levels, but the simplified design procedure differs only slightly 
from the more complicated: procedures. 


Conclusions 


1 The resistance of an absorption liner to a sound wave is 
composed of a configuration resistance and a nonlinear re- 
sistance dependent on sound pressure level. At levels higher 
than about 100 db, the nonlinear resistance controls the ab- 
sorption characteristics of the liner. This nonlinear resistance 
can be correlated with the sound pressure level of the waves 
being absorbed. 

2 The effect of increased sound pressure level on absorp- 
tion depends on the physical dimensions of the resonator. 
For most practical configurations, an increase in the sound 
pressure level will cause an increase in the theoretical absorp- 
tion coefficient to a value of 100 per cent between 120 and 150 
db followed by a continuous decrease. 

3 The effective length of the air mass vibrating in the 
resonator aperture is dependent on the per cent open area of 
the absorber configuration. Absorption liners having open 
areas ranging between | and 12 per cent should satisfactorily 
eliminate screeching combustion under most conditions. 


Nomenclature 

a = distance between the centers of individual resonators 
necks arranged in a square array, ft 

c = velocity of sound, fps 

d = aperture diameter, ft (unless noted otherwise) 7 

db = sound pressure level, 0.0002 dynes/em? = 

f = frequency, sec™ 

fo = resonant frequency, 

g = gravitational constant, 32.2 fps? 

= imaginary quantity, 


lets = effective length of aperture mass, ft 


1.0 
t 
| 
| | 


LL = absorber backing distance, ft (unless noted otherwise) - -X ~~ == specific acoustic reactance ratio for resonator array, di- 
p = excess pressure in sound waye, lb/ft? 7 mensionless 
P; = sound pressure amplitude of incident wave, lb/ft? 
Q = measure of the absorption bandwidth (Eq. [14]), dimen- Subscripts 
sionless c = critical conditions 
s = particle displacement, ft max = maximum value 
t = plate thickness, ft (unless noted otherwise) min = minimum value 
u = particle velocity, fps 
a = absorption coefficient, dimensionless 
a = absorption coefficient at resonance, dimensionless References 
ag = absorption coefficient with non-normal incidence, dimen- f , 
sionless 1 Ingard, U., ‘‘On the Theory and Design of Acoustic Resonators,” .’. 
. Z Acoustical Soc. Amer., vol. 25, no. 6, Nov. 1953, pp. 1037-1061. 
6 sa effective length correction factor, ft (unless noted other- 2 Zwikker, C. and Kosten, C. W., ‘‘Sound Absorbing Materials,’ Els 
wise) vier Publishing Co., N. Y., 1949, chap. 7. 
Ant = nonlinear correction factor, ft 3  Beranek, L. L., ‘‘Acoustics,’” McGraw-Hill Book Co., Inc., N. Y., 195), 
6 ad specific acoustic resistance ratio for resonator array, di- 4 London, A., ‘The Determination of Reverberant Sound Absorptic:: 
4 mensionless Coefficients from Acoustic Impedance Measurements,”’ J. Acoustical S: 
6) = specific acoustic resistance ratio at resonance, dimen- Amer., vol. 22, no. 2, March 1950, pp. 263-269. 
sionless 5 Ingard, U. and Labate, S., ‘‘Acoustic Circulation Effects and t! 
Nonlinear Impedance of Orifices,’’ J. Acoustical Soc. Amer., vol. 22, no. 


= wave length, ft March 1950, pp. 211-218. 
6 McAuliffe, C. E., ‘‘The Influence of High Speed Air Flow on the B 
havior of Acoustical Elements,’’ M.Sc. Thesis, Massachusetts Institute 
Technology, 1950. 
7 Morse, P. M., ‘‘Vibration and Sound,’’ McGraw-Hill Book Co., Inc , 
N. Y., 1948, second ed., p. 242. 


viscosity, lb/sec-ft 
p = density, lb-sec?/ft‘ 
o 
¢ 


= fraction of open area of resonator array, (2d?) /(4a?) 
= angle of incidence of non-normal waves 


ri 


Technical Notes 


At a symmetric body stagnation point, the application of 
the well-known similarity transformation 
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i Specific Heat Inequality Effect in the 
Chemically Frozen Stagnation Point 


_ Boundary Layer 


GEORGE R. INGER! 


(k = 0 two dimensional, k = 1. axially symmetric) 


University of Michigan, Ann Arbor, Mich. " , 
niversity o ichigar ich p d) 


Pe 
T HAS often been assumed, for purposes of mathematical U = Uf'(n) ; ' 
simplicity in the dissociated laminar boundary layer prob- 
lem, that the specific heats of the atomic and molecular species pt = Pulle 
are equal.? This, of course, is not necessarily true and may 
lead to a significant error in the temperature gradient and 
heat transfer calculation for a highly cooled dissociated bound- 
ary layer. It is the purpose of this note to give a solution to 


reduces the laminar boundary layer equations for a dissociated 
air mixture to the following set of ordinary differential equa- 
tions in the absence of homogeneous chemical reactions 


the energy equation for chemically frozen flow at the stagna- ai S.fer’ + zp” = 0 [1] 
tion point which takes this specific heat inequality into ac- 
count and explicitly shows the resulting effect on the heat + + + = 0 [2] 
conduction and total heat transfer to a highly cooled wall for 
arbitrary Lewis number. ante 
F = chemically frozen flow 

S. = Schmidt number « 

May L, = Lewis number (S, = Pr/L,) 

? The inequality of atomic and molecular species specific heats, boundary layer 

the explicit effect of this specific not Cy = Com aaze(Cy, Cym), frozen specific heat of thi 
in the results given for temperature distribution or heat transfer. mixture 


Eprtor’s Note: The Technical Notes and Technical Comments sections of ARS JouRNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are usually published without editoria! 
review within a few months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page 
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It has been assumed in Equation [2] that the atomic and 
molecular specific heats, Cp4 and C>,,, respectively, are un- 
equal constants. The boundary conditions for these two 
equations in the absence of wall mass transfer are: zp(©) = 
6p(o) = 1, & = T/T, given, and either z7(0) = 0 (com- 
pletely catalytic wall) or zp’(0) = 0 (completely noncatalytic 
wall). 

Under the assumed conditions of a highly cooled wall 
(A, <1), it has been shown by Lees (2)* that it is sufficiently 
accurate to take f(n) as the Blasius function for the purpose 
of solving the composition and energy equations. Conse- 
quently, these equations become linear and can be readily 
integrated. Equation [2] can be solved as a first-order equa- 
tion in zr’ to yield 


= 0.47 S,1/8 f, exp (-s. f, fdn )dn os 
© ? 


2Fnoncat. 1 


With the composition distribution known, Equation [2] can 
now be solved to indicate the effect of the specific heat in- 
equality term on the frozen boundary layer temperature dis- 
tribution. For the case_of zero atom recombination at the 
wall surface (zg = 1, Cp’ = 0), we see immediately that 
Ci4 — Cpm has absolutely no effect on the temperature dis- 
tribution. However, in the opposite extreme of complete 
atom recombination at the wall [z-(0) = 1, Cp’ ¥ 0], the 
specific heat difference introduces a dependence of the tem- 
perature distribution on the atomic species distribution. As 
a first-order equation in 6’, Equation [2] can be integrated 
twice to the following solution 


Ov at. = Ow + (1 — Oy) ) [4] 


where 


) —(1+Le) 
= Ba exp (—Pr san) an [5] 
P 


Now for a catalytic wall 


ip ia) 


where 
~ 1.49(0.47 S,1/3)(0.47PR) [8] 


Equation [6] clearly shows that the effect of Cp, — C>,, in- 
creases the catalytic wall temperature gradient and raises the 
entire temperature profile above the corresponding noncata- 
lytic wall solution. For the typical example S, = 0.50, 
Pr = 0.70(L, = 1.40), aa, = 0.50 and C,,/C>, = 1.40, 
Equations [6 and 8] give = Com 1.357, 
i.e., a 36 per cent increase in the wall temperature gradient 
owing to the composition influence introduced into the energy 
equation when Cp, = 1.4Cpn. 

The total heat transfer rate to the wall Q, can be written in 
the following nondimensional form in the similarity plane for 
chemically frozen flow 


prebtre(1 k)(dU, / AX) stag. 
C,(0) 6p'(0) + L,zpr’(0) 
Com ComTe 


where 


atoms 


is the average heat of formation (dissociation energy) in the 
dissociated external flow. In the noncatalytic wall case 
(zp’ = 0), there is no diffusion and only the heat conduction 
transfer exists 

GP = OATPR' + — = [10] 
In this case, the specific heat inequality does not change the 
temperature gradient but does increase the thermal conduc- 
tivity for constant Pr, hence increasing the total heat transfer 
proportionally to aa(Cp4/Cpm — 1). On the other hand, 
the catalytic wall accepts the full heat of recombination by 
diffusion [zr(0) = 0, zr’(0) = 0.47 S,'/*]; the thermal con- 
ductivity at constant Pr is unaffected by Co, — Con[Cp(0) 
= Cpym], but the wall temperature gradient is increased ac- 
cording to Equation [6]. Consequently, we have a 


Il — 


h 


varies from zero at the wall to a value of 0.20 at the edge of 
the boundary layer for the typical case a4, = 0.50, Cp4/Cpm 
= 1.40. Hence it appears permissible for most practical 
problems to use the approximation 


Cy(n) 
bal ~1- (1+ L,aszr 1) 


Concerning ourselves only with the effect on wall temperature 
gradient, then, this simplifying approximation gives 6,’(0) 
from Equations [4 and 5] as 

0.47 Pr3(1 — 


O) ~ 047P + [6] 


’ Numbers in parentheses indicate References at end of paper. 
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conduction diffusion 


In addition to the diffusion term, the catalytic wall heat 
transfer is also larger than that for the noncatalytic case be- 
cause of the relatively greater effect of the specific heat 
inequality on the conduction term. For example, in the 
previously cited case S, = 0.50, Pr = 0.70, aa, = 0.50 and — 


Cpa/Cpm = 1.4, we have 


(Gr) Cpa = Com |noneat. conduction 1.200 


(Gr) Cpa =Cpm conduction —~ 1 


The magnitude of these ratios suggests that the specific 
heat inequality effect on the heat transfer in the highly cooled 
dissociated boundary layer on a hypersonic body should be 
taken into account in engineering calculations. 
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Transient Heat Conduction in a Melting 
Finite Slab: An Exact Solution 


R. W. SANDERS! 
Lockheed Missiles and Space Div., Sunnyvale, Calif. 


Transient heat conduction in a melting finite slab is 
studied. Solutions are determined for which the melting 
rate is inversely proportional to the amount of unmelted 
material. The solutions obtained demonstrate that a 
better measure of accuracy for any numerical procedure 
is its ability to predict the melting rate and not the grosser 
quantity—the total amount melted. 


N GENERAL, the only exact solutions for the transient 

heat conduction equation are for linear problems on regions 
of simple shapes. When nonlinear problems are considered, 
solutions for arbitrary heat flux input can be obtained only by 
numerical methods. Since very little is known about the 
convergence and stability of numerical procedures for non- 
linear problems, it is desirable to have exact solutions of 
typical problems for comparison. This paper presents such 
a solution for a nonlinear transient heat conduction problem 
of current interest. This solution will be used to compare 
the effects of two different heat flux inputs on the quantities of 
interest. 

The temperature distribution and the melting rate of a 
one-dimensional finite slab, with constant thermal properties, 
subjected to a specific heat flux input will be determined. 
The slab, originally extending from z = 0 to x = a, has the 
face x = a insulated and has heat flowing in through the other 
face at a rate H(t) per unit area. Only the phase of the 
problem that involves the actual melting will be considered. 
It is assumed that all molten material is removed immedi- 
ately upon formation, so that the face which was initially at 
x = 0 moves inward, and at time ¢ is at the position x = 
S(t) (1).2. In terms of the specific heat c, density p, thermal 
conductivity K, latent heat of a fusion L, melting tempera- 
ture 7’, and the temperature distribution at the start of 
melting [7'o(z) |, the process is described by 


oT oT 

k— x <a, 
cp <2 <a,t>0 [la] 

T(x,0) = T(x) O<2 <a [1b] 
oT = 0 z=a,t>0 [le] 
=T, t>0 [ld] 
Hy(t) = - = S(),t>0 


The essential difficulty is that the problem is defined for a 
region having an unknown moving boundary. Following the 
guide of Landau (1), the moving boundary can be eliminated 
by the transformation § = (a — zx)/[a — S(t)]. If, in addi- 
tion, the dependent variables are nondimensionalized by 
defining 


t) = (T(E, t) — 


the problem can be stated as 


1 — sds 


s(t) = S()/a = [2] 
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= H(t) O<E<1,t=0 [3b | 
00/dE = 0 é=0,1>0 [3c | 

a(1,t)}=0 t>0 [3d | 


dt 


r = K/pca* [4 | 
For an arbitrary H(¢) and a nontrivial initial temperature dis- 
tribution, this problem must be solved numerically. Since: 
we are attempting to develop a solution that may be used t» 
check the accuracy of numerical techniques, we need no' 
study the direct problem but will consider the indirect on: 
That is, if ds/dt is chosen, then Equations [8a, 3b, 3c and 3: 
form a well-posed boundary value problem with Equatio 
[3e] as the specifying equation for H(t). A similar procedui 
was followed by Baer and Ambrosio (2) for the semi-infinit 
slab. If s(t) is chosen such that, with A as an arbitrary res! 
constant 


(1 — s)?= 1 — 4A*rt [5] 
the bracketed term in Equation [8a] is a constant, name!’ 
24°. After changing independent variables by setting y - 
Equation [3a] becomes 


2 2, 


Since the problem has been defined for positive t, the fact tha! 
ds/dt is nonzero at t = 0 presents no difficulties. 
The problem can now be solved by separation of variables. 


Let = ‘Then 


where \ is an unknown constant. Making the change of 
variable y? = Z changes the equation for f into a well-known 
form having two independent solutions in terms of thie 
confluent hypergeometric function 


fly) = BiF\(—); 4; y2) + — IS] 


The second term of Equation [8] can be eliminated by apply- 
ing boundary condition [8c]. Use of boundary condition 
|3d] gives the eigenvalue equation 


Fy(—2; $; A?) = 0 91 


Here 


The function ,F;(—\; 3; A%) behaves like cos (A V4 + 1) 
for large values of A. Hence, there exist \ = X, such that 
Equation [9] is satisfied. From the form of the ,F; function, 
we can show that A < 0 cannot satisfy Equation [9]. That 
\, must be real will be shown in the discussion of the orthog- 
onality property of these functions. 

To study orthogonality, change 
from fn(y) to gn(y) as follows 


only) = ~ 3; y?) [10] 
Then ¢,(y) satisfies the Sturm-Liouville equation 
+ =0 [11] 


In this form it is easily seen that the eigenvalues are real, anc 


the dependent. variable 


; that the ¢, are orthogonal over the interval (0, A). That is 
(0,m#n 12] 


A 


(Ni2,m=n 
and the normalizing factor NV, can be calculated. 
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From Equations [5 and 7] we see that g,(é) ~ (1 — 4A 2rt)>*. 
The indirect problem is linear; we can therefore write the 
general solution by superposition 


6(y,t) = — y2) [13] 


The values of the D, can be determined by applying the initial 
condition and using the orthogonality property of the ¢n(y). 
It is apparent from the relation of the D,, to 6(y) that a proper 
choice of @(y) will greatly simplify the numerical calculations 
required. With the D, determined, Equation [8e] specifies 
H,\t), the heat flux required to generate the melting rate 
indicated by Equation [5]. 

When the function (—1),Fi(—o; 3; A%E?), where Xo is the 
least eigenvalue, is chosen for the initial temperature dis- 
tribution, the general solution given by Equation [13] re- 
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Fig. 1 Temperature profiles in a melting slab 
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Fig. 2 Heat flux and melting rate vs. time 


duces t 


= 2doiFi(1 — Xo; 3; A2)- —— (1 — 4A + 


Several 
for A: 
A 
0.5 
0.85 
9.99 


» a one-term solution. This function chosen for 
6(¢) is well-behaved for almost all A. The one-term solu- 
tion for the temperature distribution, with the required heat 
flux input, is 


t) = —(1 — 4A*rt)Fi(—Ao; 3; A%E*) [14a] 
Takk 


a 
2LK A? 


- (1 — [14b] 


ca 


solutions have been obtained using various values 


Xo — Xo; 33 A) 
1.0 2.0 

403 0.5 2.3491 
0.40514 2.4933 
0.39922 2.5041 


Aspects of the last two solutions indicated are presented 
in the figures for typical thermal properties. Fig. 1 shows 
the temperature distribution at times ¢ = 0 and 20 sec. The 
initial temperature distribution is essentially the same for 
A = 1 and V0.99. Fig. 2 shows H(t) and ds/dt vs. t for 


A=1. 


Fig. 3 shows the percentage differences of H(t) and 


ds/dt for the two solutions as functions of time. These 


results 1 


ndicate that a perturbation of H(t) is reflected im- 


mediately in ds/dt even though it might go unnoticed in 
s(t). Thus a better measure of accuracy of any numerical 


procedur 
quantity 
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of the Speed of Propagation of Turbu- 

lent Combustion 


R. E. BOLZ' and H. BURLAGE Jr.2 


ve ‘ase Institute of Technology, Cleveland, Ohio 


HE AUTHORS wish to discuss Professor Shchelkin’s 

recent criticism (1)* of the combustion model and method 
of data processing employed in their paper published in 1955 
(2). The experiments discussed in that paper involved the 
ignition, by means of a single spark pulse, of a jet of homo- 
geneously mixed natural gas (methane) and air. This 
produced a small spherical flame kernel which was carried 
along by the moving stream. From successive photographs 
of the growing kernel, taken with a high speed camera and a 
shadowgraph optical system, the apparent flame propagation 
speeds were determined by direct measurement. The actual 
flame speeds relative to the unburnt gases were calculated 
by eliminating the growth rate of the flame kernel due to the 
constant pressure adiabatic increase of volume. This was 
done using 


Sactual = NRT 


number of moles, gas constant and absolute 
temperature of unburned gases, respectively 

b burned gases 

S flame speed 


where 


N,R,T = 


By calculating the constant C in this manner it is implied 
that the flame front or combustion zone width is thin com- 
pared with the radius of the flame kernels, and, therefore, 
the combustion is complete within the confines of the ob- 
served flame front. Professor Shchelkin, however, dis- 
agrees with this assumption and believes that although the 
assumption is valid for laminar flow, it is not valid for turbu- 
lent flow. He believes that for the turbulent case the com- 
bustion zone will be of the order of 5 mm in thickness, a value 
quite large compared with the kernel diameters that were 
measured in these experiments and which never exceeded 
30 mm. As a result, he predicts that the values of the tur- 
bulent flame speeds given in (2) are only about half of what 
he estimates they should be. His estimate is based upon a 
combustion situation where v’ the average turbulent pertur- 
bational component of flow velocity is very large compared 
with S, the normal flame speed. 

Professor Shchelkin comments further that the dependence 
of the flame velocity on the radius of the flame kernel, shown 
in Figs. 7, 14, 15 and 16 of (2), are not to be expected. He 
reasons that due to the damping of the turbulence with dis- 
tance one should actually expect a decrease in flame speed, 
and that this expectation would be realized if the value of 
the constant C were evaluated as he suggests, i.e., the value of 
C calculated on the basis of a thick flame of the order of 
5 mm. 

The authors seriously disagree with Professor Shchelkin’s 
criticisms as presented in (1), and this disagreement is based 
upon several arguments which will be discussed here. 

In the combustion experiments under consideration, a 
spark pulse ignited a small but finite mass or kernel of homo- 
geneously mixed fuel and air. This ignition process took 
place in a matter of microseconds, and the initial shape of 
the flame front was essentially independent of the mild 
stream turbulence as is shown in (2). The combustion sur- 


Received May 19, 1960. 

1 Professor of Aeronautical Engineering. 

? Associate Professor of Aeronautical Engineering. 

3’ Numbers in parentheses indicate References at end of paper. 


1032 


face thus be gun was dlaiiaas as the ren front moved int 
the unburned mixture and as the high temperature gas 
produced within the kernel expanded. Since there is 1 
dependence upon turbulent or diffusion mixing of the fuel an 
air (which were premixed), it appears reasonable to hy 
pothesize that for this case of very mild turbulence, where t! 
value of the root-mean-square perturbation velocity w 
usually less than 3 fps, the combustion will propagate at t! 
normal laminar flame velocity at each point on the flan 
front. Consequently the mild turbulence that was generat: 
in the flow field by screens would have only the effect 
slowly distorting and wrinkling the essentially spherical fla: 
front. An important fact here is that in the experiments 
(2) the complete time lapse for each sequence of phot 
graphs was usually less than 0.01 sec, the same order 
magnitude as the characteristic time of turbulence, and durin 
this complete time lapse, the average perturbation velocit 
component was only about 2 fps. The effect of such turb) 
lence, then, cannot be immediate, but must develop slow! 
during the very short time period of observation. Furthe: 
more, the authors can visualize no mechanism for this mi! 
turbulence to produce the thick flame front pictured |) 
Professor Shchelkin. It appears much more reasonable 
hypothesize that, in the case of the experiments of Bo 
and Burlage, since the normal flame speed is approximate! 
equal to the mean turbulent perturbational velocity compo- 
nents, the flame should progress as stated here, as an esse! 
tially laminar flame front which is being slowly wrinkle 
by the turbulence. This wrinkling is a transient proces: 
it is building up during the short time of the photographi 
observations. If the studies could extend over a longer tim 
period than that available to the experimenters (and witho: 
decay of turbulence), then one would expect an asymptoti 
steady-state value of flame propagation rate to be reache: 
This is much the same as the mechanism pictured by Scurloc 
and Grover in (3). With turbulence decay, one would ex 
pect an increase to some maximum value of apparent flam 
propagation rate due to increased flame front surface are: 
(but not a change in the basic mechanism of propagation 
followed by a decrease in this rate which would eventuall: 
approach the laminar value. 

Furthermore, based upon this model of a wrinkled lamina 
flame and the theory of (3), an equation relating the kern« 
radius (reduced by the factor C) with the time was derived 
in recent work by the authors‘ taking into account the 
damping of the turbulence. The theoretical values obtained 
from this equation were compared with new and much mor 
precise and extensive data taken in the past few years by 
Bolz and Burlage [similar to that in Figs. 7, 14, 15 and 1! 
of (2)], and the comparison showed a remarkable correlatio1 
in all cases between the theory and the experimental dat: 
as reduced by the method of the Bolz and Burlage paper 
(2). Such a correlation does not prove the validity of either 
the theory or the data reduction method but does lend sup- 
port to the authors’ arguments. 

The authors will certainly agree with Professor Shchelkir 
that if the rms turbulent perturbational velocities v’ wer 
very large compared with S, the normal flame velocit) 
(v'>>S,), the rapid oscillation of the mixture at the flam 
front would almost certainly result in a thick flame zone 
This is the case pictured by Professor Shchelkin and th 
assumption upon which he bases his flame front thicknes 
calculation in the absence of any other information. How 
ever, it is just this fact, that the v’ and S, involved in the 
work of (2) are of approximately equal magnitude, that th 
authors believe substantiates their view. 

As a further substantiation for the authors’ methods 


~ 4This work will be published shortly as an NASA publicatior 
and may be obtained on loan in prepublication form from Pro 


fessors Bolz and Burlage, Case Institute of Technology. 
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